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| We describe a universal element P in the group algebra of symmetric groups, whose 

characters provides the counting of quarter and eighth BPS states at weak coupling in 
J\f = 4 SYM, refined according to representations of the global symmetry group. A related 
projector V acting on the Hilbert space of the free theory is used to construct the matrix 
O \ of two-point functions of the states annihilated by the one-loop dilatation operator, at 
finite iV or in the large iV limit. The matrix is given simply in terms of Clebsch-Gordan 
coefficients of symmetric groups and dimensions of U (N) representations. It is expected, 
by non-renormalization theorems, to contain observables at strong coupling. Using the 
stringy exclusion principle, we interpret a class of its eigenvalues and eigenvectors in terms 
of giant gravitons. We also give a formula for the action of the one-loop dilatation operator 
on the orthogonal basis of the free theory, which is manifestly covariant under the global 
symmetry. 
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1 Introduction 

The half-BPS sector of M = 4 U(N) SYM has been a very fruitful area of study in 
the context of AdS / CFT [H [2J, [3] from a number of different perspectives. It has shown 
how different classes of half-BPS local operators in the large A" theory can be mapped to 
different types of physical objects in space-time : Kaluza-Klein gravitons, strings, giant 
gravitons and classical space-time geometries [21 HI El El El El El CD]- The study of the 
two-point function has been central. It defines an inner product on the space of states 
corresponding to the local operators. The diagonalization of the inner product in the 
space of half-BPS operators in terms of Young diagrams [7] has provided the tools for 
identifying the physical objects in different semi-classical limits. For example we expect a 
state describing two distinct giant gravitons to be orthogonal to a state of a single giant 

Ed- 

Recent years have seen progress on the quarter and one-eighth BPS sectors. From 
the space-time point of view, there is a better understanding of the geometry of giant 
gravitons and their connections to simple harmonic oscillator systems [TTJ [121 CHS US] • 
From the point of view of gauge- invariant operators, we have a diagonalization at finite 
N for free Yang Mills i.e with the g\ M = HQ EZl HB1 US] ■ It is known that there is 
a jump in the spectrum of BPS operators from zero to weak coupling. It is conjectured 
that there is no further jump as the coupling is tuned to the strong coupling region [T3] . 
This makes it extremely important to understand the weak coupling sector in complete 
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detail. The lessons we learn from these sectors will be useful for an understanding of the 
sixteenth BPS sector where we can access black holes with finite horizon areas. 

Henceforth when we refer to BPS operators we have in mind quarter BPS operators or 
bosonic eighth BPS operators (more precisely eighth-BPS operators with bosonic lowest 
weight states) at weak coupling. The quarter BPS problem involves two matrices X±,X2 
transforming in the adjoint of U(N) and a U(2) global R-symmetry which rotates these 
two matrices. The eighth BPS problem involves X±, X 2 , X 3 and a U(3) global R-symmetry 
which rotates these. In our considerations below, we can keep the global symmetry general 
in the form U(M). The general U(M) case does not have a direct application in M = 4 
SYM, but may have applications elsewhere in theories with M copies of adjoints. 

The explicit construction of BPS operators was studied in [201 121] • A class of quarter 
BPS operators involving a small number of X\,X2 were constructed in the SU(N) the- 
ory at weak coupling. In [TBI [T7], the diagonalization of the two-point function of BPS 
operators in the free (zero-coupling) U(N) theory at general finite N was accomplished. 
The solution for the case of a general number n of X±,X2 was given in terms of the 
group theory of S n . A Fourier basis constructed using Clebsch-Gordan coefficients and 
matrix elements of S n was found to accomplish the diagonalization. The one-loop BPS 
operators can be characterized as being orthogonal, in the free field theory metric, to 
the descendants [21]. Based on [22] it was conjectured that these BPS operators have 
non-renormalized two- and three- point functions. The idea of characterizing BPS op- 
erators using orthogonality was developed in a finite N set-up, using the zero coupling 
inner product to give expressions for the BPS operators [16]. The expressions were based 
on the notion of dual bases of operators [231 I2H 125] and arrived at expressions involving 
inverse dimensions of U(N) representations, which will play a role later in this paper. 
A missing ingredient was an explicit description of the descendant operators. Exploiting 
some developments in BMN operators [26J a characterization of the BPS operators in the 
1/N expansion was given in [2Z] using the inner product in the planar limit. For related 
recent developments on BMN operators, see [28J. 

In this paper we describe the symmetrization operation on traces as a linear operator, 
which we call V, acting on the orthogonal Fourier basis for the free theory. This allows 
us to characterize the descendants as the kernel of V. We use the symmetrization matrix 
to derive a formula for the two-point function of the BPS states to all orders in the 
1/N expansion. A manifestly finite N construction is also given. It relies on an infinite 
dimensional Hilbert space which contains the Hilbert spaces for any N. The geometry of 
intersections of subspaces in this Hilbert space plays a crucial role. 

The paper is organised as follows. We will review the relevant background in Section [21 
It will become clear that an explicit description of the linear operator V implementing 
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the symmetrization will be important. The operation of symmetrizing traces appears in 
the simpler problem of counting BPS operators. In Section [3] we review the standard 
counting results and consider a refined counting according to representations A of U(M). 
This is done with a new approach based on a universal element in C(S n ) related to trace 
symmetrization, which we denote as P, and whose characters provide this counting for 
any U(M). We derive a generating function which gives the structure of P. This turns 
out to have a rich combinatoric structure, which we describe in Appendix [El We show 
in Section [4] that analysing the fine structure of P leads naturally to the symmetrization 
operator V. This allows the construction of the BPS operators in Section 14.21 and leads 
to the matrix of two-point functions in Section 14.41 Section gives a description of 
the finite A" construction. This discussion relies on a careful description of an infinite 
dimensional Hilbert space Ji spanned by states |A, Ma, R, t) containing the labels of the 
free orthogonal basis. An inner product which we call the inner product plays a key 
role. The finite A" Hilbert space can be realized as a quotient of % by the image of a 
finite A" projector. The Soo inner product induces a simple inner product on the finite 
N Hilbert space, which is consistent with the finie A" cutoffs. Clarifying the interplay 
between the finite A" projector and V, viewed as an operator in % is the key to the 
finite A" construction. In Section El we analyze the properties of the matrix of BPS two- 
point functions. Using the stringy exclusion principle [2HJ EI] , we interpret some of its 
eigenvectors in terms of giant gravitons and their spacetime excitations. 

We collect miscellaneous technical material in the Appendices. Appendix [B] also gives 
a new i7(3) covariant formula for the mixing of descendant operators under the action 
of the one-loop dilatation operator, which has a similar structure in terms of symmetric 
group Clebsch-Gordans as the symmetrization operator V . 

2 Review and Remarks on BPS operators 

We review here recent work on the construction of quarter and eighth (bosonic) BPS 
operators. It is known that the spectrum of BPS operators at zero coupling M = 4 SYM 
(the free theory) is different from the spectrum at weak coupling. At zero coupling, any 
gauge invariant operator constructed from holomorphic combinations of the three adjoint 
complex scalar fields Xi, X 2 , X 3 is BPS. At weak coupling only those operators in the 
kernel of one- loop dilatation operator %2 are BPS [32J. In [T6l [T7] , the diagonalization of 
the BPS operators in the free theory at finite A" was accomplished. We will begin with 
a short summary of this work, where a Fourier basis was found to be central to the free 
field diagonalization. We will follow it with a summary of observations from [2U [111 [27J 
on the construction of BPS operators at one loop. This will motivate the construction 
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of an operator V on the Fourier basis, which implements symmetrization of traces. In 
Section [3] we will obtain new results on the counting of BPS operators, which provide an 
independent avenue and concrete steps towards the operator V . 



2.1 The orthogonal basis in the free theory 

In diagonalizing the free two-point function [161 EL it was extremely useful to exploit 
the notion of Fourier transformation as applied in the context of the symmetric group. 
To set the stage, let us recall that the Fourier transform relation between a circle S l and 
the set of momenta Z can be understood in group theoretic terms as a transformation 
relating the group U (1) to the space of inequivalent representations Z. The transformation 
from permutations of n elements to representations is a generalization of this Fourier 
transformation to the symmetric group S n . 

Gauge invariant multi-matrix operators can be written in the form 

_ tr n (Xg a) _ h in 

In the first expression, the operator = X ai (g) X a2 • • • <g> X an where the matrix X ak 
can be viewed as a linear map from the iV-dimensional vector space V to itself, with 
matrix elements (X afc )^ fe . The index a transforms in the fundamental M-dimensional 
representation of U(M). By tr n we mean a trace of operators in V® n , i.e tr n (A) = 
^il'-'-T- The permutation a takes integers l---n to ot(l) ■ • • a{n). In the context of 
representation theory, the a label for operators naturally Fourier-transforms using matrix 
elements Df-{a) to a set of labels R,i,j, where R runs over irreducible representations 
(irreps.) of S n and i,j run over states in the irrep. The state indices a of the U(M) 
symmetry can be transformed using Schur-Weyl duality into labels A, Ma, m where A 
labels simultaneously irreps of U(M) and S n , Ma labels a state in the irrep A of U(M), 
and m labels a state in the irrep A of S n . The transformation is achieved by the Clebsch- 
Gordan coefficients C^ M k for the decomposition 

v m = Q v u(M)^ v s n (22) 

A 

This leads us to consider 

OA,M A ,m,R,i,j = ^ A^( a )^A,M A ,fc^g,a ( 2 - 3 ) 
a,a 

However the labels a, a are a redundant description of the gauge invariant operators. 
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Imposing the correct invariances leads to operators 



Oa,Ma,R,t - -^T^ E S ^ f™ T D ij( a )^A,M A ,mPs,a (2-4) 



a, a 



where S ^ f^ 1 T are Clebsch-Gordan coefficients for coupling R R <S> A to the one- 
dimensional irrep. of S n . These operators also have the virtue of diagonalizing the free 
two-point function. We will call this the Fourier basis of the free theory. The inverse map 
giving the trace basis in terms of the Fourier basis is 

Oa, a = E V^RC^ lA ' m D^a)S R J^O AMA ^ T (2.5) 

A,Ma,R,t 

We defined c^' MA ' m = {C^ M )* which satisfy standard orthogonality relations for 
Clebsch-Gordan coefficients [T7] . 
Defining 

qA,m a ,r,t _ rr~ r ,A ' M ^' m n R (^\Q R R A ' r 

^a,a - y (1 R L 'd L> ij\ OL ) b i j m 

7T (2-6) 



^a'm a ,R,t — j ^A,M As m^ij( a )^ i j m 

we may write a more compact version of the transformations (12. 4p (12. 5p between trace 
and Fourier basis 



We have 



and 



Ua,a - S q UK,R,M h ,r 
^A,R,M A ,r — ' A,M a ,R,t 

^A 2 ,R 2 ,T2,M' A 



(2.7) 



y^^AiJV/A.A.r^q.q ^ ~ ^^M^, M'^ S Rlj R 2 5 T1 jT2 (2.8) 

E ^ AAT -^2r AAT = E ^ WM/rVW) (2.9) 

A,M A ,/?,r 7 

The 5 at is defined as 

djv (<*) = >. ; 2.10 

R: Cl (R)<N 

When n < N, this is the sum over all irreps of S n , which is equal to 



5(a) = 1 for a = identity element of S„ 
= otherwise 



(2.11) 
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The right hand side of (12. 9p is the identity operator in the finite N Hilbert space, which 
we will denote as 1^' and return to in Section [5j 

The two-point function defines an inner product on the gauge-invariant operators 



<Oi|C? 2 > = lim a:f((O 1 ) t (x 1 )(O 2 )(0)) (2.12) 

XI— >oo 



Here and throughout the paper we use the zero- coupling but finite N two-point function 
and the associated inner product, unless specified otherwise. The inner product on the 
trace basis evaluated by Wick contractions is: 



7,0-G.S'n 

76SV 



(2.13) 



On the Fourier basis it is diagonal [T7] 



nlDimi?! 



xm(n) 



(2.14) 



—, ^Ri,R2^t 1 ,t 2 5a u a 2 5m a ,m' 

"Hi 



There is no d\ as in eq. (110) [TF] because we defined our Clebsch-Gordans to be the 
coupling from R £g> R £g> A to identity rather than R £g> R — > A. In addition because of 
the \fdn in the definition of the operator of d R in the denominator rather than d 2 R . We 
explain the derivation of the result (12. 14ft in Appendix [Fl 
The Q factor 

VL = ^N c °~ n o- (2.15) 

is a central element in the group algebra C(S n ), i.e it commutes with any element of S n . 
Defining the operator J 7 on the Fourier basis as 

, a 2 ,m' A2 ,r 2 ,t2 nlDinxRi 

IAuAIa^RuT! - °Ai ,A 2 °M Al ,M'a 2 °n ,r 2 Cfli ,R 2 , 

-xx xx XrA* l ) 

— " Ai ,A 2 "Ma, ,M' a t1iT2 RiiR2 — 

2 



We see that 



(Oa u m Ai:RuTi \0 A2Ma R2tT2 ) - {F) k 2 \ m K a '!r 2 \t 2 ( 2 - 17 ) 
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The operator Q is defined with inverse matrix elements 

^)a 2 ,m' A2 ,r 2 ,t 2 - d M m hi Al ,M' A2 d T1 ,r 2 d Rl , R2 - } — — 



' (n-i\ ( 2 - 18 ) 
x xx XrA* 1 ) 

0\ 1 ,A 2 0M Al ,M A2 <Vi , T2 dijj ,R 2 — 



The second line is only meaningful in the region n < N, where Q can be inverted. 
The equation for the inverse dimension in terms of characters XRl j^ — - is useful in two- 
dimensional Yang Mills theory (see e.g. equation (6.6) in [53]). In the leading large N 
(planar) limit, where n <C N, the two-point function approaches the identity 

(Oa 2 ,Ma 2 > r * ' t 2 I > M Ai < R i > r i ) ~~ ^ ,M'a 2 < r z > r 2 I ®M,M Al ,Rl >Ti ) planar = #Ai ,A 2 ^M Al ,Ma 2 ,r 2 Sr 1 ,R 2 

(2.19) 

and we may write 

(Oa 2 ,M'a 2 ,R2,t 2 \Oa 1 ,Ma 1 ,Rut 1 ) = (OA 2 ,M' A2 ,R 2 ,T 2 \J I 'OA l! MA 1 ,Ri,T 1 ) P l&na,r (2.20) 

Since the large N limit of Q is 1 the 2-point function (12.131) in the trace basis approaches 

(^|a>)planar = £ ^/(^ W) ■ (2.21) 

yes„ 

This shows that the 2-point function in this limit is zero unless the two operators have 
the same fields in the same trace structures. This is large N factorization. 



2.2 Towards the construction of BPS operators 

In this section we review earlier work on quarter and eighth-BPS operators at one-loop 

[221 EH HH [27]. 

These operators, which have vanishing one-loop anomalous dimensions, are annihilated 
by the one-loop dilatation operator % 2 [311 132] : 

■H 2 O bps = 0. (2.22) 

where I-L2 in the U(3) sector is given by 

Tfa = -i teK , (2.23) 

Equivalently, the protected operators are precisely those which are orthogonal, in the 
inner product defined by the zero-coupling two-point function, to the SUSY-descendant 
operators C D : 

(0 D |0 BPS ) = 0. (2.24) 
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Indeed, from hermiticity of I-L2, we have 

= (O any \n 2 BPS ) = (H 2 C any |C BPS ) = (C D |C BPS ) (2.25) 

It is expected that the 2- and 3-point functions of these operators receive no further 
corrections at higher loops. 

The descendant operators are easily found: they are those which, when written as a 
gauge invariant product of single traces, contain a commutator of fields [2T] : 

O d = tr([X ai ,X a2 ]X a3 X a4 . . .)tr(. ..)... (2.26) 

This is true even when the total number of X's denoted by n is in the region n > N, 
in which case, since there is no unique way to write an operator as a product of traces, 
an operator would be a descendant if there exists a way to write it with a commutator. 
Since we have an explicit set of O d , the task of finding O b is just finding the orthogonal 
subspace according to f)2.24p . This has already been pointed out and used in [21] . The 
procedure, however, turns out to be complicated because of the non-planar inner product 
and does not provide a convenient way to express (9 BPS when n grows large. 

In order to make progress [271 126] it is useful to use (12.201) in order to express the 
non-planar two-point function as 

(0»\0 BPS ) = (0 D |^ BPS ) planar , (2.27) 

Now pick the space of operators O s which are orthogonal to O d in the planar inner 
product 

(O D |O S )planar = 0. (2.28) 

The space of protected operators can be written as 

e> BPS = go s . (2.29) 

where Q is defined in (12.181) as the inverse of J 7 . It follows simply from (12.271) that such 
operators will be orthogonal to O d in the non-planar two-point function 

(0 D \gO S ) = (0V£O S >planar = (0° 1 S >planar = (2.30) 

and thus they are protected. The action of J 7 , which acts diagonally by multiplication 
with "j^ 1 ^ in the Fourier basis is the operation 

T tr(X s a) = tr(X s Q a) (2.31) 
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where Q is the central element in the algebra of C(S n ) defined in (I2.15p . Indeed, using 

_ Xnjty rr r A,MA,m n B { w R R A;t n 

~ d R 5 U ij\ a > b ijm U A,M A ,R,r 

yA,M A ,m nR/n \qRRA;t, 



= VdRCr iA ' m D«(na)S" t ?^0 A , MA , R , T 
= O s , Q a (2.32) 

To go from the second line to the third we use the fact that Q is a central element, hence 
proportional to the identity in any irrep. R by Schur's Lemma. Similar steps can be done 
with the Q matrix. In the large N limit, where n < N, we can define as an element 
in C(S n ) and 

Q tr(X 3 a) = tr(X s IT 1 a) (2.33) 

We now want to make a final point, which was not fully expressed or exploited in the 
earlier literature. The space of operators O which are orthogonal to O d in the planar 
inner product is the space spanned by products of symmetrized traces. The symmetriza- 
tion operation replaces each single trace with the sum over all permutation of elements 
inside the trace 



symm [tr(X ai X a2 . . . X a J ] = — ^ tr(X aa(1) X a<r(2) . . . X a ^ n) ) = Str(X ai X a2 . . . X a J 

(2.34) 



n! 

U&Sn 



The operators O s are linear combinations of 

Str (X ai X a2 . . . X ak ) Str (X bl X b2 ...X h ).... (2.35) 

The orthogonality to descendants can be seen, using the property of large N factorization, 
as follows. An operator O , orthogonal to descendants, is a sum of individual multitrace 
operators. Take one multitrace term of it 

tr (X a X b A) B, (2.36) 

where A is a product of matrices and B is some remaining multitrace operator. Now the 
whole O s has to be orthogonal to a descendant 

tT([X a ,X b ]A)B. (2.37) 

But in order for that to happen the O s must also contain a term 

tr (X b X a A) B, (2.38) 
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so that 

(tr {{X a X b - X b X a )A) B | tr {{X a X b + X b X a )A) £?) planar = 0. (2.39) 

Of course, it can be that tr (X b X a A) is already the same as tr (X b X a A) either because 
of a = b or because of cyclicity of the trace, but that does not change the fact if one is 
contained in O s then the other must be as well with the same coefficient. Therefore, since 
every two adjacent matrices in O have to be symmetrized, we conclude that the whole 
O must be built from fully symmetrized traces as in (12.351) This allows us to conclude 
that the space of one-loop BPS operators is spanned by 

O bps = G [Str (X ai X a2 ...X ak ) Str (X bl X b2 ...X bl )...]. (2.40) 

In this paper we will develop (12.401) and arrive at a description of the BPS operators 
in the 1/N expansion in terms of an explicitly defined matrix QV in the free Hilbert 
space of SYM. The matrix V, related to symmetrization, will be described in Section I4~2"l 
Incidentally, the properties of V will lead to a one-line proof of (12.281) in Section 14.31 We 
will derive the two-point function on the BPS operators as a matrix {VQV). In Section [5] 
we will extend these results to finite N, where we will identify a projector Vi, which 
approaches V in the large N limit. In the next section, we will derive some new results 
on the counting of BPS operators, which will lead naturally to the operator V. 

3 Counting and the universal element in C(S n ) 
3.1 Chiral ring counting 

Here we briefly review some known results about counting of BPS operators. 

The expected number of bosonic gauge invariant eighth-BPS operators with given 
R-charges (ni, 712,713) can be calculated using either the Chiral ring counting arguments 
[T3] or alternatively using the description of the BPS sector as a Hilbert space of N 
3D harmonic oscillators [33]. In the chiral ring gauge invariant operators made from 
Xi,X 2 ,X 3 are counted modulo equivalences following from 

[X 1 , X 2 ] = [X 2 , X 3 ] = [X 3 , Xi] = 0, (3.1) 

It is argued that this leads to a counting of traces of diagonal matrices [13] : 

tr ((X i r(X 2 ) b ^X 3 r) tr ({X 1 ) a >{X 2 ) b >{X 3 )<») . . . tr ((X 1 )°*(X 2 ) 6 *(X 3 ) Cfc ) (3-2) 

or as a symmetric polynomial of the N eigenvalues of the three matrices. This counting 
of states also arises from the collective coordinate quantization of the coordinates of the 
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BPS moduli space determined by (13.11) upon reduction to + 1 dimensions on S 3 [35]. In 
the limit of large N the counting is given by the partition function Z(x\, x 2 , x 3 ): 

00 ^ 00 

Z(xi,X 2 ,X 3 ) = ■ —— y, Z niri2n3 xl X X2 2 X3 3 , (3.3) 

i+j+fcX) 12 3 rii,n2,rt3=0 

with numbers Z ninan3 providing the number of operators with charges (711,712,713). This 
will be valid as long as rt\ + n 2 + 713 < N. 

For a £7 (iV) gauge group at finite N, the counting is given by the generating function 
Z N (xi,x 2 ,x 3 ) which can be expressed as: 

Z(v;x 1 ,x 2 ,x 3 ) = - - = 1 + y^Z N (x 1 ,x 2 ,x 3 )v N (3.4) 

ij,fc=o ^1^3X3 JV=1 



OO 



Zjtffei, X 2 , X 3 ) — Z^. niTl2na x" 1 X2 2 ^3 3 (3-5) 

ni ,712,713=0 

In terms of counting this is equivalent to limiting the number of traces k appearing 
in the operators such as (13.21) to at most N (this is certainly not true in the actual 
construction of BPS states at finite N as can be easily checked using the examples of 
Appendix O). This limit arises from the fact that there are only N distinct eigenvalues 
to build the symmetric polynomials, or, in the 3D harmonic oscillator picture, from the 
fact that the collection of labels (a*, bi, c,) indicate the excitation numbers of k particles, 
and we are limited to at most N excited particles. 

For the purposes of this paper we would like one extra refinement, that is, to distinguish 
operators in different representations A of the global symmetry group £7(3), which rotates 
the three matrices X a . We will represent A = [Ai, A2, A3] by a Young diagram with at most 
three rows of lengths Ai > A2 > A3 of U(3). Operators transforming in the representation 
A contain a total number of matrices equal to the number of boxes in the Young diagram. 

Ai + A 2 + A 3 = m + n 2 + n 3 = n (3.6) 

We denote by the multiplicity of A. Each representation A contributes to the full 
partition function the Schur polynomial in variables Xi, therefore the partition function 
can be expanded as [361 EZ] : 

Z(xi,x 2 , x 3 ) = 22 -^a Xa(xi, x 2 , x 3 ). (3.7) 

A 

This can be inverted to calculate Ai\ (see for example (A. 23) in [37]): 



M 



Z(x x ,x 2 ,x 3 ) J 



Xj Xj t 
i<3 



(3.8) 



r *l+2 A 2 + l A3 
L l ^2 3 
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where the square brackets with the subscript is an instruction to take the coefficient of 
the indicated monomial from the function inside. Equivalently we can define 



1 ~ x\ ) ( 1 ~ xW ( 1 ~~ ~£~ 2 



(3.9) 



and write 



•^[Ai.Aa.Aa] 



M(xx,x 2 ,x 3 ) 



(3.10) 



Organizing the operators into representations of U(3) works equally well for finite N. We 
can define M. n-,a as the multiplicity of the A representation in the space of operators at 
finite N. We can again define 



M N (xx,x 2 ,x 3 ) = Z N (xx,x 2 ,x 3 ) ^1 - ^1 - ^ ) 
and write for finite N multiplicities A^at ; [a 1 ,A2,a 3 ] 



X3 

x 2 



(3.11) 



M 



W;[Ai,A 2 ,A 3 ] 



M N (xi,x 2 ,x 3 ) 



\*2 



(3.12) 



We refer the reader to (ID.2[) in Appendix [D] for some examples of 

We can also refine the counting for specified U(l) charges ni,n 2 , n 3 into those belonging 
to specified U(3) representations A. Let Z£ in2n3 , Z^. nin2Tl3 be these refined multiplicities 
at large and finite N respectively. They are given by 



ninans = M ^ #( A , [™l] , M , M) 

M N . A g(A, [nx], [n 2 ], [n 3 ]) 



N\n\n-2.nz 



(3.13) 



where g([ni], [n 2 ], [n 3 ]; A) standard group theoretic numbers, namely the Littlewood- Richardson 
coefficient for combining three single- row Young diagrams of lengths nx,n 2 , n 3 to give the 
Young diagram A (see for example [37]). 

The representation multiplicities M.n-,k thus provide a refinement of the counting 
7 7 

Zv nin2ri3; Zj 'N ' \n\n2nj, • 



7 -^T 7 A 



«2»3 



J N;n\n2n^ 



Z~2 Z N;n! 



(3.14) 



«2«3 



In the rest of the paper we will focus on these representation multiplicities .Ma, M- 



iV;A- 
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3.1.1 Universality from counting formulae 



The above arguments can be run for the case of quarter BPS operators with U(2) sym- 
metry. We have the partition functions 



Z(xi,x 2 
Z N (x 1 ,x 2 ) ■- 



OO j 

n \ _ x i x i = x ™ 1 

i+j>0 1 2 m,n 2 =0 



™«2 
x 2 



oo 

TT — 

. . „ 1 — VX\ In 



" ni ,712=0 



(3.15) 



x 
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and representation multiplicities 

-M[ Al ,A 2 ] = z i%i, X2){xi - x 2 ) 
Mn;[\ 1 M]= Z N (x 1 ,X 2 )(x 1 - X 2 ) 

These can be usefully expressed in the form 



Ai+l A 2 
r l x 2 



(3.16) 



- A/i [A 1 ,A 2 ] 
■MjV;[Ai,A 2 ] 



^Ai,A 



^Ai+1,A 2 -1 



^JV;Ai,A 2 ~~ ^ 



(3.17) 



iV;Ai+l,A 2 -l 



An important point is that these expressions in fact follow from the ones for £7(3) simply 
by setting A3 = 0. If we consider generating functions with M variables, as appropriate 
for U (M) symmetry, for sufficiently large M, they contain all the information about the 
multiplicities Ai\ for any Young diagrams A having n boxes, with n < M. Of course 
when M < n, we are only interested in Young diagrams with less than or equal to n rows. 
All Young diagrams A with n boxes correspond to the complete set of representations 
of S n , so we may expect structures entirely in S n which contain information about the 
multiplicities Ma, -Mn-,a for all A with n boxes which are irreps. of S n . In the next 
section, we will show that A^a is related to a universal element P and M.n-,k to a finite 
N version thereof pW. 



3.2 Universal element 

We will now present an alternative method of counting the BPS operators and multiplic- 
ities A^a using the symmetric group techniques. We will develop a notion of universal 
element ¥ in the symmetric group algebra whose characters give the multiplicities. The 
element is easy to define but the precise coefficients of different permutations in P (which 
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only depend on the conjugacy class of these permutations) are not easy to infer directly 
from the definition. We will find that the generating function of BPS states Ai(x\, X2, £3) 
in ( 13.1 Op can, after a change of variables, be used to derive a generating function for these 
coefficients. Analysing the structure of P will lead to the construction of BPS operators 
in Section |U 

First consider the space of abstract strings of n operators X\, X2, X3 without putting 
any traces: 

X ai ®X a2 <g> ...®X an , (3.18) 

where a\ take values 1 to 3. Most of what we will discuss is true for any global symmetry 
group U (M) with indices running over 1 to M, in the case at hand we will specialize to 
M = 2 for the quarter BPS sector and M = 3 for the eighth-BPS sector. Each "letter" X a 
can be thought of as a state in the vector space Vm, which is the fundamental of U(M). 
The whole string is then a state in the tensor product representation V^ n . 

There is the natural action of the symmetric group S n on these states where the 
permutation a shuffles the indices: 

^D^M^.^I^. (3.19) 

By Schur-Weyl duality the tensor product representation can be decomposed into irre- 
ducible representations of both Z7 (3) and S n : 

V® n = © A Vf (M) ® Vfr . (3.20) 

This immediately tells us that if our operators were defined by ordered strings of X a 
as in ( I3.18P then the multiplicity of the U (M) representation A would be just d\ - the 
dimension of V A n . 

In our case of interest the operators are instead products of traces like (13. 2p . where 
the order of operators inside the trace does not matter. Let us still view them as states 
in V^ n , but with extra structure, for example: 

tv(X ai X a2 )tv(X a3 X aA )tv(X a5 ). (3.21) 

The groups U(3) and S n act on the indices just like before. The only difference now 
between counting states like (I3.2ip and those in ( I3.18P is that if we change the ordering 
of the operators within a trace or if we swap all operators between two traces of the same 
length it is still considered to be the same state. This restriction can be formulated in 
terms of the S n action: in our example we want to consider the states which are invariant 
under permutations a = (12), a = (34) and a = (13) (24). 

Consider a general trace structure which we represent by a partition p = [l Pl , 2 P2 , . . .]. 
This means there are Pi traces of length i, and we have £V ipi = n. We will first count how 



16 



many states there are with the given p and later sum over it. The set of permutations 
which leave the state invariant given the trace structure p will form a subgroup of S n 
which we denote G(p). It consists of all permutations within traces: 

GM = XiiSt)*, (3.22) 

with each S* C S n understood to be embedded in S n in such a way to permute its 
corresponding trace. Regarding S n as the permutations of integers 1 ■ ■ • n, the 1-cycles 
correspond to S'f 1 acting on 1 • • - pi, the 2-cycles correspond to Sf 2 acting onpi + 1 • • -pi + 
2p2- G(p) also contains permutations which exchange traces: 

G 2 (p) = x^, (3.23) 

where each S Pi permutes all the elements in the pi traces of the same length %. The full 
group which leaves the state invariant is the one generated by G\{p) and G^p), which 
can be seen as a semi-direct product^]: 

G{p) = G t {p) x G 2 (p). (3.24) 

For completeness note that: 

\Gi(p)\ = Um p \ |G 2 (p)|=n^) ! ' \G(p)\ = \G 1 (p)\\G i (p)\. (3.25) 

i i 

Note that we have chosen a specific embedding of G(p). Different permutations a in the 
conjugacy class p give different embeddings G(a) of the same group in S n . We can now 
define a projector which will project onto the subspace of V^ n invariant under G(p): 

; 6G(P) (3.26) 
= |Gi(p)||G 2 (p)| ^ ^ ° T ' 

Choosing different a also leads to different ¥ a , a special case of which is P p : we will return 
to this freedom in the next section. The element P p is exactly what we need to account 
for the multiplicity of distinct operators with the trace structure p: the dimension of the 
space projected to is just the trace of the projection operator 

M p = tr v «»(Pp). (3.27) 



3 It is a semi-direct product in the sense that G2 provides an automorphism of G\ by conjugation in 
the S n - giGxg^ 1 = Gi for 92 £ Gi. Then group multiplication of two elements in G\ x Gi can be defined 
as (51,32) • (hi,h 2 ) = (gi(g2hig2 1 ),g2h2) where gi,hi e G% and 32,^.2 G C? 2 - 
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What we in fact want to know is how many times a particular U(M) representation A 
appears in this space. For that we also project onto the and divide by its size: 

Ma ' p= Dn^V tr ^"(W' ( 3 - 28 ) 

where DimA is the dimension of U(M) representation A. Because of the Schur-Weyl 
duality f)3.20p we know that Pa just projects onto the M ) (g) Vf n factor, and so: 

Ma ' p= Dm^V (M) ^ (Pp) 

= ffiA tr ^^^ (3 ' 29) 
= Xa(P p )- 

Here xa is the symmetric group character and this is an expression entirely in terms of 
S n quantities! The second line follows from the fact that P p consists of permutations, 
therefore, it acts only in the V^ n factor. Then the trace factors into U (M) and S n parts, 
with the trace over identity in the U(M) canceling DimA. This result relates nicely to the 
fact we mentioned earlier, that for states living in V® n , the multiplicity of representation 
A would be just = XaOO- Instead here because of the invariance of products of traces 
under some permutations captured by P p we lose some states, which is reflected in counting 
by taking trace of the projector P p instead of identity. 

We are just a short step away now from the formula for the total multiplicity .Ma- All 
we need to do is to sum over all partitions p describing different trace structures. Also 
note that the precise form of P p depends on how the G(p) group is embedded in S n , i.e. 
how the traces are ordered. However, the counting does not depend on that, and so we 
want to make things more uniform by summing over all embeddings. That gives us the 
final form for what we will call the universal element P in the algebra of the symmetric 
group C(5 n ): 



P 



£ (h £ 




(3.30) 



n\ \G(p)\ 

Note that P depends on n but we will leave that dependence implicit. 
Now the formula for counting multiplicities is just 



M A = xa(P). (3.31) 
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Just to emphasize, this is exactly the same .Ma as in ( I3.8p . but now we have expressed 
it as a character of an element P of C(S n ). The first nice thing about this result is 
that the same P will give multiplicity for any U(M) representation A. When the global 
symmetry is U(2), we are only interested in A = [Ai, A2], i.e in Young diagrams with at 
most two rows; for U(3) we apply the formula for A = [Ai,A2,A3]. Second, the P itself 
is an interesting object to study. It is universal in the sense that it is defined uniquely 
for every S n and there is no additional data that goes into it. Since P contains averaging 
over conjugation by 7, its terms will be sums over conjugacy classes, and we will find that 
the coefficients weighting these conjugacy classes are integers (after multiplication by n\) 
with interesting properties. 

Before we go on let us provide a short example of P and .Ma's to get a better feeling 
for them. Take n — 4, states made from four matrices. Using (I3.30P we can explicitly 
calculate: 

P(S 4 ) = ^ (3S [4] + 3E [3)1] + 7E [2)2] + 7E [2il)1] + 15£ [liliU] ) , (3.32) 
where S p is the sum over all elements in the conjugacy class p: 

S p = ^(7. (3.33) 

Now if we go and compute the characters of P we get: 

A* [4] =X[4](P) = 5 
%=X[3,i](P) = 2 (3.34) 

M[2,2]=X[2,2](P) = 2. 

These are indeed the right multiplicities. Converting to the more familiar numbers of 
operators for different (n 1? n 2 , n 3 ) according to (13. 7p gives Z i00 = 5, Z 3W = 7, Z220 = 9. 

Note that xa(Pp) has some vanishing properties which are easy to understand by 
thinking about symmetries and trace structures. When p is a single cycle, then the 
character vanishes for any A with more than one row. When it has two cycles, the 
character is zero for any A with more than two rows. This reflects the well-known fact 
that the simplest quarter (eighth) BPS operator has two (three) traces at large N. 

We refer the reader to Appendix [D] for explicit examples of P at various n. 

3.3 Generating function for P 

As already mentioned, P will be an important ingredient in the construction of operators 
in this paper, so we would like to know more about it than the formula (I3.30p . from which 



19 



it is hard to see what the sums evaluate to. First, we find that it can always be written 
as 

nS n ) = -J2t P Z p , (3.35) 

71. 

phn 

where the coefficients t p are integers. These are numbers assigned to each possible parti- 
tion p, and, therefore, have various combinatorial properties, which might be of interest 
in themselves. For certain simple sequences of partitions they belong to integer sequences 
tabulated in [38J . We discuss them further in Appendices [Dj [E] and provide more explicit 
examples there. 

Consider an exponential generating function t(y) for t p in the sense that: 

Ey Pl y P2 y P3 
tp=[iPi2P23P3...] J -r^-r^T ■ ■ ■ (3.36) 

so that power of yt indicates the number of cycles of length i in the partition. We claim 
that this generating function is in closed form: 

t(y) = exp ( jr^. ( e d ZiV« - l) J = JJexp Q ( e d ^^ - l)\ (3.37) 

Vd=i / d=i 



or more explicitly 

t{y) = exp ^( e «i+W+W+... _ l) + I ^+2y 4+ 2y 6+ ... - 1) + I ( e 3y3+3 y6+ 3, 9+ ... \ 

(3.38) 

This is as close as we can get to an explicit formula for P. 

We will prove that (13.371) generates the coefficients of P correctly by verifying the 
counting formula (I3.3ip with the assumed t(y) against the known result f)3.10p . The 
expression xaQP) can be evaluated directly from the generating function t(y) with the 
help of Frobenius character formula. It states that a symmetric group character can be 
written as 1371: 



Xa(p) 



H( Xi - x 3 ) H (xf + x\ + . . . + x k n ) Pk 



(3.39) 



n — 1 + n — 2+A2 
c l x 2 
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where p is a partition. Using this we can express xaQP) as: 



phn 



phn 



t. 



(3.40) 



Q(x; - Xj) 



Li<j 



|Sym(p)| fe 



Q(x* + x^ + ... + x*) 



x l x 2 



In the second line we used the fact that E p contains n!/|Sym(p)| elements all in conjugacy 
class p, where Sym(p) is the automorphism group of permutations in p and 



|Sym(p)| = n^ ! >' Pi - 



(3.41) 



Let us now perform a change of variables 

1 



We get 



y k = -{x\ + x k 2 + ...+ <) 



i<j 



phn \Sym(p)\ k 



\k Pk {y k 



n-l+Ai ra-2+Ao 



PA' 



,*<J 



phn k 

~\i x i - x i) x t(y 



n — 1 + Ai n — 2+Ao 



_i a ."-l+^l x "-2+A 2 ^ 



(3.42) 



(3.43) 



with exactly the generating function t(y), with the understanding that the y k variables 
have to be substituted according to (I3.42p . 

We can evaluate the function t(y) in terms of x k variables: 



V i=l/ L d=l L V k i 



[X: 



(3.44) 



Let us simplify this using Y2 k =i( xk )/k = 1°§(V(1 — x )) anc ^ some other infinite series 
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relationships: 



i (f) = exp|gi^exp^log T -^- 



exp 



oo 



. d=l 



l x l x 2 x 3 • • V 

i+j+k+...>0 



(3.45) 



E log ( 

j+k+ ...>o Vl-(xlx 

^— 



,-y» J ,-yt't' "\ 

x 2 x 3 . . .J 



i+j+/c+...>0 



We find, in fact, that t(x) in terms of Xi variables is precisely the partition function 
Z(xi, x 2 ,x 3 ) from (13.31) only with arbitrary number of variables 



t(xi, x 2 , x 3 ; x i>3 = 0) = Z(xi,x 2 , x 3 ) 



(3.46) 



Evaluating (I3.43P for the specific case of counting U(3) representations A = [Ai, A 2 , A 3 ] we 
then immediately find 



X[Ai,A 2 ,A 3 ]( 



Z(x 1 ,x 2 ,x 3 ) Y\_( x i - x j) 



i<j 



(3.47) 



A-,+2 A 9 + l i, 
j, x j, 2 ^3 



This is precisely A^[a 1 ,a 2 ,a 3 ] according to (13.81) ! That means we have explicitly shown 
Xa(P) = -Ma and, therefore, proved (I3.37p . 

3.4 Finite N counting 

We can also introduce a finite- N version of the universal element, denoted as F^ N >, and 
defined as 



p(w) = E £ t^t" 1 ) 



(3.48) 



The superscript on the summation symbol indicates that the summation over par- 
titions p = [l Pl 2 P2 ---] is restricted by Y2Pi — N. This constraint follows from the 
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description of the chiral ring counting in terms of N particles in a simple harmonic oscil- 
lator, which can be described as the counting of products of symmetrized traces with no 
more than N factors. By repeating the arguments explained in deriving (I3.3ip . we obtain 

-Mat ; a = Xa(P W ). (3.49) 
As in (13.351) there is an expansion in conjugacy classes 

P(A ° = ^[E^ (3.50) 

Let us introduce conjugate variable v to N again and define a deformed generating func- 
tion t(u; y): 




(3.51) 



Then the picking off the coefficient of v N in t(y\ y) will provide a generating function 
tjq{y) for the coefficients t p in P^: 

t(iy;y) = J2^t N (y). (3.52) 

JV 

This can be proved by noting first that the derivation (I3.45P with t(u;y) instead of t(y) 
gives M.{y \ x±,x 2 , x 3 ). Then we see that 

XA(F (iV) ) = [*a(PM)U 

= [t(v, y)(xi - x 2 )(x 1 - x 3 )(x 2 - X3)} uNx x 1+ 2 x x 2 +i x x 3 (3.53) 

= M N -A, 

where by P(z/) we just mean the object generated by t{y;y). This is again in agreement 
with Chiral ring counting results. 



4 From counting to operators on the Hilbert space 

The counting of BPS operators in Section |3] used a universal element P G C(5* n ) defined 
in (13 .30 p . It is a sum of projectors P p over conjugacy classes p. Equivalently it can 
be viewed as a sum over the symmetric group of projectors P a , where a G S n is a 
permutation specifying a concrete embedding of the cycle structure p. We examine this 
in Section 14.11 and also derive a related operator p. The operator p is then used in 
Section l4~2l to construct an operator V which acts on the free basis Oa,m a ,R,t and produces 
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linear combinations of symmetrized traces. With the help of V we then construct the BPS 
basis valid at N > n. We then pause to discuss the properties of V in Section 14.31 and 
conclude with the calculation of the two-point function on the BPS states in Section 14.41 
In Section we will also see how the operator V allows the construction of BPS 
operators at finite N. 



4.1 A p- map from £(S n ) -> C(S r 

As we remarked in Section [3] the element F p in C(S n ), which counts the number of U(M) 
representations A among trace structures specified by the cycle structure p, depends on a 
choice of embedding of the group G(p) in S n . Different permutations a in the conjugacy 
class p are associated with different groups G(a). We can define an element 

P « = |G^)| £ " <"> 
There is a related linear map p : C(S n ) — > C(S n ) 

p(a)=F a = J2Pd, a P (4-2) 

P 

A trace tr n (X^a) is symmetrized by the action of F a 

symm[tr n (X 3 a)] = ^ p /3>Q tr n (X /3(S )a) (4.3) 

P 

We extend the definition to sums of traces by linearity, i.e the symmetrization of a sum 
of multi-traces is the sum of terms obtained by symmetrization each multi-trace. 
The universal element in the previous section f)3.30p is a linear combination 



P 



n. 

phn 7SS- 



2^ ~\ Fq (4-4) 



n 



\- |Sym(q)| 

- — ~\ — 

We can define the coefficients in the Fourier transformed basis for C(S n ). 

P?g^^^(^)p^A» (4-5) 

a,P 
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R2 




Rz 




R 2 



Figure 1: Diagram for p 



Then 



P/3,. 



R,S,i,j,k,l 



dR<k 
n\ n\ 

dp d. 



Dim Psl\D^a) 



(4.6) 



P(«)= £ ^>^A»^(/3)/3 (4-7) 

R,S,i,j,k,l 

Figure [3] gives a diagrammatic expression of the equation (14.51) . On the right hand side 
we use the standard diagrammatic form for operators. On the left we we have short-hand 
which we will use subsequently. 



4.2 The V operator and the BPS basis 

In this section we proceed with the explicit construction of the BPS operators, annihilated 
by the one loop dilatation operator "H 2 - As explained in Section |2T2| in order to do that, 
we have to consider operators of the form (I2.40p . The construction involves two tasks: 
symmetrization and the application of Q. A priori it is not clear how to accomplish both in 
a way to arrive at a closed form expression. The reason is that if we start off with a naive 
basis of multitrace operators, the symmetrization is easy to implement, but the action of 
Q (I2.33P is complicated, involving the inverse Q^ 1 which is not well defined for n > N. 
On the other hand if we start with the diagonal basis of the free theory 0a,m a ,r,t, the 
Q action is straightforward (I2.18p . but it is not clear how to impose the symmetrization. 
Here we will show how the p^^ defined in the previous section allows us to do just that. 
The procedure described in this section is illustrated in the Appendix [C] with concrete 
examples of A = [2, 2] and A = [3, 2]. 

We define the following basis for the BPS operators 



(4J 
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In order to evaluate this expression, the key step is to calculate how the symmetrization 
acts on 0^m a ,r,t- We use (14. 3 p and the fact that it acts linearly: 



symm [O a>Ma ,r, t ] = ^ f™ ^ ( a )^A,M A ,m s y mm t tr ( X a a)] 



-Y. S * fm T D«(a)Ci MAtm pp, a tr(Mt) «) (4.9) 



n! 



R 



ni 



S , 1fi' r A?(a)^(/3)Cl )MA)m ,p /3 , a tr(X g a) 



Now re-express the trace in terms of the Fourier basis, use orthogonality of the U (M) x S n 
Clebsch, the fusion property of the S n Clebsch : 



Symm [0 AjMa ,R,t] = ^ 5 1 f m T D ii( a ) D mm'{P) C tM K ,m>Vp, 

fl A 1 ,M Al ,mi n R 1( ^ Q Ri Ri A i;n 



En „Ai,iM Al ,mi n R 1( „\ Q Ri Ri Ai;n n 
V «Ri G a U ab \ a ) b a b mi Ai,M A ,.Ri,"rL 

Ai,A/ A ,fii,ri 



E V d R d Ri \^ nR R A, r r)R( n \ rjA ffl\nflif„\ n q Ri Ri A;n n 
1] » J m L) ii\ 0L ) L> m rax\P) L> ab (aJP/3,c^ a 6 mi U A,Mk,Ri,ti 

Ri,ti a,ft,S,Ts 

E V d R d Ri \^ qRR^rqRS Ri,t sq RS Ri,rs n S( n \ n A ( o\_ 
-j 2^ d S b i 3 m b ika b jlb U kl\ a ) U mmi\.P)Vp,a 

Ri,t\ a,/3,S,T S 

q Ri Ri A;ri ,-/-» 
X D a mi ^A.Ma^i.ti 

E Vd R d Rl d S ^ s k j q R r a,t q R S Rur S Q R S Ri,t s qRi Ri A,n n 
, PAmmi J i j m ° i fc a ° j I 6 ° a 6 mi ^A,A/ Al iii,ri • 

S,T3,Rl,Tl 

(4.10) 

In the final line we recognized the sum over a, /3 as a Fourier transform of the p^ 
This allows us to define a matrix P: 

A2,Ma 2 ,R-2,T2 

\l JAi,M Ai ,Bi,ti 

_ r r d S \Jd Rl d R2 Rl g R2 . T g kl k2 Rl g R2 . T Rl Rl A . T1 R2 R2 A . T2 

— 0Ai,A 2 0M Al ,M A2 n \ °h ki ji P A m A m A D i 2 fc 2 h D '1 »2 m A °jl j 2 m A 

SjT m A ,m A 

(4.11) 

and write the symmetrization as 

Oi,M A ,R,r = symm [O a>Ma , r>t ] = (^)£K; T1 <?A,M A ,Ki,n- (4-12) 

-Rl,Tl 
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We have thus been able to express the symmetrization of individual traces as the trans- 
formation of the free orthogonal basis by a matrix V. 

Let us make a remark about the notation. Since V is diagonal in the Ma label and 
also independent of it, we will often use a shorthand matrix (P)^ 'jji'n understood as V 
without the Sm Ai ,m a factor: 

On the other hand, even though V is diagonal A, it does depend on it, so we will keep the 
A label. The same applies to the matrices T and Q. Note also that whenever a repeated 
A or Ma index appears, the summation convention does not apply. 

Evaluation of (14. 8 p is now straightforward. Using ( I4.12p the BPS operators can be 
written as 

°KLat = E (^tt^CW^n (4.14) 

fll,TL 

or more explicitly 

The equation (I4.14p is a key equation for BPS operators and we will devote the rest 
of the paper to analyzing these BPS operators in detail. We have written an expression 
for the space of operators annihilated by the one-loop dilatation operator 

%05L,H,r = (4.16) 

in terms of known group theory quantities. They are nicely packaged in V ( 14. lip involving 
Clebsch-Gordan coefficients and pp )Ct , and in the Fourier basis of operators (12. 4p for the 
free theory. To be more precise, the above space of operators provides an overcomplete 
basis. The labels of {A, Ma, R, t} uniquely label the free basis, but we know that the 
kernel of "H2 has lower dimensionality. A linearly independent basis can be found by 
analysing the inner product on the the BPS operators. We will give an expression for the 
inner product in Section 14.41 and discuss its eigenvalues and eigenvectors in the context of 
the dual space-time physics in Section [6j 

The expressions (14.141) . (I4.15P should be understood as valid for n < N . They already 
contains hints on the problems we might encounter at n > N. The Dimi? factor in the 
denominator will be whenever Ci(R) > N making the basis ill-defined. This is related 
to the fact that Si -1 appearing in (12.331) is not well defined as n > N. We will solve the 
finite problem in Section [51 
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4.3 Properties of V 

The matrix V will be a crucial object throughout the rest of the paper and so we pause 
here to note a few of its properties. In addition to this section we provide further analysis 
of V in Appendix |A] 

• The operator V can be seen essentially as the operator pp, a transformed appropri- 
ately into the basis of free operators labeled by |A, M\, R, r). As is apparent from 
the definition ( 14. lip , the transformation is done by applying four Clebsch-Gordan 
coefficients. We can also visualize this transformation by expressing it diagram- 
matically (see Figure [2]). We can also write the matrix V as a sum over conjugacy 
classes [a] in C(S n ), which replaces the sum over R in (I4.12p . See Appendix lAl 

• The matrix (V)^'^ 1 '^'^ is proportional to 5a 1 ,a 2 5a/ Ai ,a/^ and is symmetric under 
exchange of R\, T\ with R 2 , r 2 . 

This means that it has the hermiticity property 

7?t = p 
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under any inner product where the (A, Ma, R, t) label an orthonormal basis. The 
free field inner product does not have this property, except in the planar limit 
(n <C N) where it simplifies to (I2.19p . In the finite N construction of Section [5] we 
will define an inner product which has this property for any n. 

Since V is constructed as the operator on the free field Hilbert space which projects 
to symmetrized traces, we expect it satisfies the projector equation 

-p 2 = -p 

We prove this equation directly in the Appendix |A] 

The hermiticity of V in the planar inner product gives an elegant way to see that 
the symmetrised traces are orthogonal to the descendants in the planar limit - a 
property we first proved in Section 12.21 States in Ker(P) are descendants. States 
in lm(V) are symmetrised traces. The hermiticity of the projector V ensures that 
Im(P) is orthogonal to Ker(P). Explicitly 

= (VO D \O s ) pl ^ r = (0 D \VO s ) pl ^ r = (0 D |O s ) planar (4.18) 

For each A, Ma, the states 0\ t M A ,R, T span a space of dimension 

Y,C(R,R,A) (4.19) 

R 

where C(R,R,A) is the multiplicity of the one- dimensional irrep. in the S n tensor 
product R x R x A. Note the symmetry under transposition of the Young diagram 

C(R,R,A) = lrY,XR{<r)XR{<r)XA{<T) 

(7 

= ^E(- 1 ) CT ^( a )(- 1 ) CT ^ T ( (T )^( (T ) (4 ' 20) 

(7 

= C(R T ,R T ,A) 

We expect that, given a choice of orthogonal basis labelled by r in Inv(i? <g> R<g> A), 
there is a corresponding choice r T in Inv(i? T ® R T (g> A), such that 

wtss = ± ( p )i;f ;| ( 4 - 21 ) 

This will be demonstrated in examples by the explicit construction in Appendix [UJ 
using bases in r space constructed as eigenstates of generalized Casimirs defined in 

EH. 
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4.3.1 Counting Check 



The counting of states in a given irrep. A of U(M) is given, using f)3.30p by 



(4.22) 




It is useful to write this as 



M A = V XK^^-APar 1 ^ 1 ) 



(4.23) 



On the other hand (P^A^n 

at fixed A, M\ is a projector matrix to the space of 
symmetrized traces. The dimension of this space should be the large- N number of BPS 
operators .Ma- In general, the dimension of the projection must be given by the trace of 
the projector matrix. We can thus perform a consistency check by calculating the trace 
of V: 



E (T) \^,Ma,R,t _ V ™ n A < R\ D s („\ a R a S Q R R S,t qR R S,t qR R A,n Q R R A,tj 

R,T T,TL,S 



Now given the form of P CT1 we know that if o<i commutes with <j\ then tr(P (7l a"2) = tr(P CTl ). 
This is then equal to M A according to (I4.23P as required. 

4.4 The two-point function for BPS operators 

Given the expression of BPS operators as (I4.14p we may easily evaluate the two-point 
function matrix. 

i,M Al ,R l! - n \ C) A2,M' A2 ,R 2 ,T2) = (GVO Klt M A ^ Rl , Tl \GVO K2tM ^ R ^ T2 ) 




(4.24) 




anar 



(4.25) 
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We have used the fact that the two-point function of Oa,m a ,r,t is diagonal in the labels 
and given by the matrix J 7 in equation (I2.17p . which is an inverse of Q. This leads to the 
formula 

(4.26) 



BPS 

M Ai ,Ri,ti 



ivgv) 



Al,i?l,Tl 
Al,i?2,T2 



for the two-point function on of BPS operators, valid at all orders in the 1/N expansion. 
This is the main result of this paper. 

Since the Clebsch-Gordan coefficients of the symmetric group can be chosen to be real 
[4*0] . the matrix VQV is a real symmetric matrix. We will develop a manifestly finite N 
construction in the next section to arrive at ViGn'Pi- m Section [6] we will use analyticity 
arguments to show that the matrix VQV can be used to extract finite TV properties. 



5 Finite N and geometry in Hilbert spaces 

The derivation of the two-point function (14.261) is valid to all orders in the 1/N expansion 
but requires N > n. Here we will describe a framework which allows the derivation of 
a formula for the 2-point function at finite N, with no restriction on n. This treatment 
contains some general insights, notably the geometry of the interplay between different 
projectors in appropriate Hilbert spaces, which should be useful in thinking more generally 
about finite iV issues, e.g in the sixteenth BPS case in M = 4 or in more general theories. 

5.1 Finite N Hilbert space as a quotient and BPS operators 
from intersections 

Start with an infinite dimensional Hilbert space H, graded by a natural number n. For 
each n, there are finitely many states. The states are labeled by |A, Ma, R, t) for each 
natural number n. A, R are Young diagrams with n boxes, Ma a state in the U (M) 
representation corresponding to A, r runs over the multiplicity C(R, R, A) which is the 
number of times the one-dimensional representation of S n appears in the tensor product 
R® R® A. We can regard the states as made of gauge invariant polynomials in variables 
(Xi)), (X2)), where i,j extend from 1 to 00. These are traces e.g tr^iA^A^). The 
relation between the traces and the group theoretic labels is given by the equation (12.71) 
involving the matrices S, T, involved in the free-field diagonalization problems solved in 
[TBI [TT] . constructed from Clebsch-Gordan coefficients of U(M) x S n , matrix elements 
and Clebsch-Gordan coefficients of S n . They make no reference to the U(N) gauge group 
hence they make sense in the N = 00 set-up we are describing here. We can define a 



31 



simple inner product on 7i 



(A', M' Ah B!, r'|A, M A , R, r) Soo = 5 a ,, a 5 Ma! m' a ,5r,r>5 t , t ,, (5.1) 

which we will call the Soo inner product. Note that its form matches the planar limit inner 
product f)2.19p . but there is a subtle difference in that (-l^s^ is defined on a manifestly 
N = oo Hilbert space H, while the planar product is defined as the limit of the finite 
large A" inner product for traces involving a total of n matrices chosen from Xi,X 2 ,X 3 
when n ^ N. 

Finite N physics can be obtained by replacing the formal variables (Af a )*- with matrix 
elements {X a ) l j of N x N matrices X a . It was shown in [36], [16] that the finite N counting 
of gauge-invariant operators can be achieved simply by cutting off R with the condition 
that its first column does not exceed N in length, i.e C\{R) < N. It is useful therefore to 
introduce a projection operator 

I^A,M A ,R,r) = h M *- R ' T) « Cl(R)£N (5.2, 
\0 if ci(R)>N 

The finite N Hilbert space is isomorphic to 'H/Kei(I^). We find that a very useful way 
to think about finite N physics is to consider projectors in %. Alongside 1^ we will also 
be interested in V, as defined by the equation (I4.1ip . In the finite A" Hilbert space, we 
do not think of V as a linear operator defined by its action on traces, because the traces 
do not form a good basis, hence the emphasis on (I4.1ip . Understanding the interplay 
between V and 1^ will be the key to understanding finite A" physics of BPS operators. 
Since V is real and symmetric in the A, R, r basis, it is a Hermitian operator with respect 
to the Soo metric. The operator is also Hermitian with respect to the same metric. 
In the following we will describe the construction of a new orthogonal projector Vi, where 
orthogonality is meant with respect to the metric. 

The metric induces a metric on the quotient H/Ker(l( N ^), which is isomorphic to 
the finite N Hilbert space / w- N>) of finite A^ matrices. On this finite A^ quotient it is 

(A l7 M Al , R\ , 7i|xW|A 2 , M A2 , R 2 , t 2 ) Soo = 5 A (5.3) 

This metric can be evaluated on the trace basis using (12. 9p to be 

7 ' (5-4) 

7 R:a(R)<N 
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For operators with n <§C A 7 ", this is just the delta function expressing large A" factorization. 
Operators only overlap if they are a product of like traces with like field content. 

Whether at large N or finite N, the descendants are in Ker(P). Whenever a single 
trace contains commutators, it is a descendant and it is annihilated by the symmetrization 
operation V. For example in the case of A = [2, 2] in the Appendix [Cj we have the 
expression 2"C[3,i] — 4|0[2,2] + |C[2,i 2 ] for the descendant in the Fourier basis. For finite 

matrices, the last term is identically zero, so we may equally write |C[3,i] — ~h%P[2,2\- 
Equivalently, without ever talking about finite A" matrices, but only about projectors in 
"H, we obtain the same expression by considering the image of -V)mU which 

contains the linear combination |C^3,i] — ^C[2,2]- We look for a hermitian (under the 
5*00 metric) projection operator Vi which will annihilate the descendants, so we have 
ViX^ n '{1 — V) = 0. We also would like this projector to map states in U {N) to H {N \ so 
that X^Vj = VilS N > = Vj. Combining these we have Wj = V{P = Vj. Together, we 
see that these equations characterise a projector for the intersection 

lm(V) D Im(Z (A °) = Im(Pf) (5.5) 

An equivalent way to arrive at this is to consider the orthogonal complement in % to 
Ker(I^) U Ker('P) in the metric. The dimension of this complement is equal to the 
number of independent BPS operators, since imposing the finite A" condition in % sets to 
zero states in Ker(/^ 7V - ) ) and setting descendants to zero imposes the vanishing of states 
in Ker('P). This orthogonal complement is Im(/^) n lm(V) so we have the orthogonal 
decomposition 

U = (Ker(/ (7V) ) U Ker(P)) © (lm(X (7V) ) n lm(V)) . (5.6) 

A natural question arises how, given two projectors and V, one explicitly con- 
structs a hermitian projector for their intersection space 

Pi = Plm(K N ))nIm(V)- ( 5 - 7 ) 

Let us consider the null space of X^ N \l — V)X^ N \ It will be composed of two orthogonal 
subspaces: the null space of X^ itself and the null space of (1 — V) inside of lm(X^)\ 

Ker(j( 7V )(l - = Ker(X^) © (im(I^) n Ker(l - V)) 

= Ker(X^) © (im(X^) n ImfP)) . 

This means if we had a projector for Ker(I^(l — V)X^) then we would just multiply 
by X^ on each side to extract the second component 

Vl = X^V KcT{xW{1 _ v)xW) X^. (5.9) 
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The question of constructing Vj is then reduced to the question of, given an arbitrary 
matrix X, constructing the projector for its null space Vkct(x)- It can be written as follows 



A- 



?W) = k - - f[(X - Aj), (5.10) 

lU=l{- A i) i= i 

where % runs over the non-zero eigenvalues of X and k = Rank(X). This, in turn, can be 
rewritten in terms of Schur polynomials 

k 

The characters x\H](X) are the Schur Polynomials for Young diagrams with j rows of 
length 1 or one column of length j. This manipulation then allows us to write the 
intersection projector Vi fairly explicitly in terms of matrix multiplications as 

k 

J ~° 5.12 
X = X^ N \l-V)X^\ 

k = Rank(X). 

Also note that using the decomposition (15. 8 j) we can express the number of operators 
Mn,a as 

M N . A = |Im(X (Ar) ) nlm(V)\ 

= |Ker(Z w (l - V)1 {N) )\ - |Ker(X (Af) )| (5.13) 
= Rank(Z^) - Rank(jW(l - V)1 {N) ). 



5.2 Finite N construction of BPS operators 

We will prove that the BPS operators at finite N are given by the image of QPj as an 
operator in %. 

We first restate the construction of the BPS operators at large N in language that will 
generalize at finite N . The large N construction starts with two key observations. The 
relation between the planar inner product (.l-)^ and the 1/N corrected inner product 
(.|.) is given by (see (12301) ): 

(Oi|O a > = (Ch\F0 2 )s„ (5.14) 
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The second observation is that descendant operators contain commutators [X ,X&] inside 
a trace. These operators are annihilated by the operator V we constructed. From these 
facts and the Hermiticity of V with respect to we easily see that the operators in the 
image of QV are orthogonal to descendants, and thus BPS: 

(o D \gvo) = (o D \Tgvo) Sac 

= (O D \VO) Sx =0 = (VO D \O) Soc =0 

To develop a finite version of this argument, we need to carefully consider the definition 
of Q at finite N . Define Qn = I^'^I^' which is manifestly zero for Ci(R) > N and well 
defined at finite N. Acting on a Hilbert space H with states |A, Ma, R, t) it is non-zero 
only on the operators obeying c\(R) < N, i.e where the Young diagram R has first column 
of length smaller than N. The second ingredient is the projector Vi defined in (I2.20p . 

Note the following useful properties of Q, J 7 , Z^ N \ V, Vi 

Qn? = FGn = X (Ar) 

(5.16) 

V{P = vv I = v I 

Consider operators QnViO. These operators are manifestly well-defined at finite N and 
have the following 2-point function with descendants 

(o D \g N VjO) = (o D \jg N VjO) Soc 

= (0 D |JWP / O) 5oo 

= {0 T° )S ~ (5,7) 
= (0 D |mO) 5oo 

= {VO^\VjO) Soo 
= 

In the first line we have used the fact that the metric of the free theory is related to the 
Soo metric by the matrix J 7 . In the fourth line, we have used the fact that V is Hermitian 
with respect to the metric, and annihilates the descendants. This shows that the 
Q^VjO construct BPS operators. 

e> BPS = v T g N o (5.i8) 



The number of these operators is |Span(Pr)|, which as explained in section [57TI correctly 
matches the finite N counting. 
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The two point function on BPS operators is the matrix ViQ^Vf. 



{GnViO^ iM ' a ,R 2 ,t 2 I GnViOa 1 ,M Al ,Ri ,n ) 
= (Oa 2i m' A2 ,r 2 ,t 2 \'PiGnJ : 'Gn'Pi Oa 1 ,m Ai ,r 1 ,t 1 )s x , 
= {VrQ^V^^ (5.19) 



M Al ,i?l,Tl 
R 2 ,T 2 

I)A 1 ,M Al ,Ri,ri 

M,R 2 ,t 2 



It is interesting to note that the geometry of intersections Im(X^) nlm(P) has played 
a key role. The projector V was constructed from the universal element P which gave the 
large N counting. Although it is easy to give indirect arguments, based on points already 
discussed, that the above finite N construction with Vj agrees with the finite N counting 
from P^), there was no direct finite N construction following directly from F^ N \ the way 
the large N construction based on V followed from P. This was not apriori obvious and 
is intriguing. 



6 The matrix ViQVi and space-time physics 

According to non-renormalization theorems [22} ETJ UT], we expect that the two-point 
functions of BPS states, constructed in this paper as the matrix ViQVi acting in the free 
Hilbert space, will be unchanged as the coupling g\ M is increased from weak to strong 
coupling. We will attempt to extract aspects of the 2-point functions at one-loop which 
are most likely to have a strong coupling interpretation. 

The set of eigenvalues, and their multiplicities, being independent of a choice of basis 
are of particular interest. In the regime n < N, the matrix ViQVi simplifies to VQV . 
We will also argue that, by analytic continuation, that information about VjQVi can 
be extracted from VQV . In particular, there is a simple way to get the characteristic 
equation of the former from the latter. The behaviour of the spectrum for n near N 
is specially interesting in view of the stringy exclusion principle (SEP) [2"9~t ED] arid the 
related property of the growth in size of giant gravitons with angular momentum [31]. The 
matching of these spacetime phenomena with gauge theory operators is well-understood 
in the half-BPS case using a Young diagram classification of the operators [7]. We will 
use the physics of the SEP to make a start towards the identification of gauge theory 
operators corresponding to BPS operators beyond the half-BPS case. We will then turn 
to a discussion of how to find a labelling analogous to the Young diagrams of the half-BPS 
case in the case of quarter and eighth-BPS operators at hand. 
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6.1 Eigenvalues of VQV and VjQVi 

The matrix VQV is derived as the matrix of two-point functions of BPS operators calcu- 
lated in the regime where N >n. As described in Section EJ when N <n the two-point 
function of the well-defined BPS operators is instead given by a modified matrix ViQ^Vi- 
However, the eigenvalues of VQV, even though calculated in the large N regime, are 
explicit functions of N. Therefore, one may ask, what happens to the eigenvalues as 
we continuously decrease N below n. For an example of the functions see ( 1C.17I) in 
Appendix O 

First, let us review the half-BPS case A = [n]. Since V = I, the BPS operators are 
just rescaled free operators 



and the two-point function 



The eigenvectors coincide then with the orthogonal half-BPS basis labelled by Young 
diagram R [TJ. The eigenvalues Xr = r ^ D ^ R are equal to N n divided by a polynomial 
[\ij(N — i + j), where the product runs over boxes in the Young diagram R, with % 
labelling the row and j labelling the column. The eigenvalue \r diverges when the first 
column has length ci(R) in the region ci(R) > N. The eigenvalue of a BPS state at 
large N, when analytically continued to a regime where that BPS state is not part of 
the finite N Hilbert space, diverges. In the original half-BPS construction or in the free 
orthogonal basis the states were normalized so that the eigenvalues continued to zero. 
In our construction of the general quarter and bosonic eighth-BPS states, the divergent 
normalization is more natural. 

Interestingly, we find the same behavior to be general for any A. That is, if we 
calculate eigenvalues at large N and then take the limit for any finite N, the divergent 
eigenvalues will correspond to the states that drop out from the Hilbert space. On the 
other hand, each finite (but non-zero) eigenvalue corresponds to a state that remains 
BPS. Equivalently, we can read off the finite N characteristic equation by factoring out a 
common denominator and dropping terms in the numerator which vanish at finite iV (see 
the discussion around flC.29j) ). In fact, the surviving finite eigenvalues of VQV with the 
corresponding eigenvectors precisely match those of ViGnVi, calculated in the manifestly 
finite N framework. In the discussion here we are not interested in zero eigenvalues of 
VQV, because they just span Ker('P) and do not depend on N. Also, since Ker('P) 
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is contained in Ker(Vi), the zero eigenvalues of VQV will always correspond to zero 
eigenvalues of ViQnVi. 

This correspondence between eigenstates of VQV and ViQnVi can be explained as 
follows. Consider the matrix VQV at a value N e — N + e approaching some finite N in 
the regime N < n. The matrix Q will contain finite and divergent parts: 

Q = Q N + 0(e~ 1 ). (6.3) 

The finite directions are those in Im(X^) and the divergent ones are in Im(l-X^), where 
the projector X^ is defined in ( 15. 2p . Now, VQV can be viewed as a metric induced by 
Q on the subspace Im('P). The divergent direction of Q will induce a divergent direc- 
tion on Im('P) which is a projection of Im(l — X^) to Im('P) and can be expressed as 
In^T-^l — X( N ')V). Therefore, as e — > there will be a set of eigenvectors of VQV with 
divergent eigenvalues which span lm{V(l —X^ N ')V). The remaining non-zero eigenvalues, 
however, should be finite and correspond to non-divergent directions. Since eigenvec- 
tors with different eigenvalues have to be orthogonal (in the metric which is used 
in the subsequent discussion), they will span a space in Im(P) which is orthogonal to 
lm(V(l — X^ N ')V) and that is precisely 

lm(V) n Im(X {JV) ) = lm(Vi). (6.4) 

That means that the finite-eigenvalue eigenvectors of VQV live in lm(Vi) and on that 
subspace VQV = ViQnVi. Therefore, the finite eigenvalues of the two operators match. 
It is also plausible, and justified by examples we have studied, that the precise forms of 
the eigenvectors of ViQnVi can be recovered as limits of those of VQV. 

As an additional evidence that interesting strong-coupling spacetime physics can be 
expected from the eigenvalues of VQV, note that the 2-point function in the half-BPS 
case is — -. The expansion of is related to conjugacy classes of S n and in turn 
to Casimirs [33j 02] [39] • In the spacetime picture they are related to conserved charges 
related to multipole moments of LLM geometries [13] . 

6.2 Identifying giant graviton states 

We now turn to a discussion of the physics of BPS states which disappear as n is reduced 
to values below N . We expect on physical grounds [31], [TTJ that there will be states in the 
region of n = N + k with k order 1 as N — > oo, which can be interpreted as single quarter 
or eighth-BPS giant gravitons along with excitations of the space-time background : single 
and multi-particle gravitons. 

In the half-BPS sector, with Young diagrams of type A = [n, 0], there are Young 
diagrams R which have a long row of order N, and a few rows with order 1 boxes. 
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These correspond to a giant graviton large in the AdS^ directions along with background 
spacetime excitations. Then there are states close to a sphere giant (large in the S 5 
directions) with one column of length order N, and a few columns with order 1 boxes. 
There are also states consisting of multiple giants, each having a finite fraction of the 
total angular momentum, along with background excitations. These correspond to Young 
diagrams with a few long columns along with shorter columns. There are also generic 
states corresponding to Young diagrams of shape far from the above types. 

More generally, we can consider a sequence of representations of the form A = [n— A, A], 
where A is order 1. We expect that a single giant graviton will have well-defined U(2) 
quantum numbers. As n increases for fixed A, the total energy of the giant will increase, 
and the energy available for excitations of the background spacetime will increase. 

The distinguishing characteristic of the states which consist of a giant expanding in the 
S 5 , along with order 1 spacetime excitations is that as k = n — N increases from 0, their 
size will grow as the angular momentum increases and they will start disappearing from 
the spectrum. In fact all the states that disappear as k increases from 0, will be precisely 
the states consisting of a sphere giant along with space-time excitations. Consider then 
the difference 

A.Mjv=n-fc;[n-A,A] = -M 7V=oo;[n-A,A] — -MjV=n-fc;[n-A,A] (6-5) 

This will measure the number of states in the irrep. [n — A, A] which disappear due to 
finite N effects. A first consequence of our discussion is that this number should become 
independent of n as n increases. We have accumulated extensive numerical evidence 
for this behaviour. We will return to these properties of &-M.N=n-k;[n-\,\] and related 
counting formulae in future work. 

Along the lines of the discussion based on continuity in section 16.14 we expect that 
states that are close to the cutoff, will have finite but large eigenvalues of ViQVi (or 
VQV). Similarly we would expect that the states with the smallest eigenvalues of ViQVi 
will be associated with single AdS giants and their excitations. 

Clearly an important open problem is to find a precise description of the eigenstates of 
ViQVi which will allow a more precise identification of the giant graviton states, notably 
to distinguish which states just correspond to single or multiple quarter-BPS giants and 
which correspond to giants along with excitations of the background spacetime. We turn 
to a discussion of avenues towards this goal. 

6.3 Orthogonal BPS basis 

Eigenvectors of the 2-point function matrix are precisely the ones that allow an identifica- 
tion to semi-classical states 0, [10] in the half-BPS case. The same line of argument has 
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been used above for states close to giant gravitons. We expect that the same can be ap- 
plied to corresponding to deformations of the bulk geometry, i.e quarter and eighth-BPS 
generalizations of LLM. 

In Section POl we have derived the two-point function on the BPS states as the matrix 
VQV: 

(osv.ft.ni °m%^r 2 ) = ^,^M Al ,Mi 2 w;;t n ( 6 - 6 ) 

Consider an eigenvector of VQV with eigenvalue A M (A) and components . The 

label v is a multiplicity label for eigenstates with a fixed eigenvalue A M . The eigenvector 
equation is : 

E (WViZ'ZK;^ 2 = KB^ lTl - (6-7) 

We know there will be Ai\;N of eigenvectors for which A^ is non-zero and finite. The 
remaining eigenvectors having zero eigenvalues are just the ones spanning the kernel of 
V and so they correspond to the descendant non-BPS operators, while the eigenvectors 
with divergent eigenvalues drop out of the finite N Hilbert space. An orthogonal basis 
for the BPS operators is then 

snBPS — rA;Riy BPS (n o\ 

R,T 

with the two-point function given by 

//nBPS I /-/-jBPS \ _ TyA;RiT 1 / T> f> / r > \A,R 1 ,T 1 R A l R 2T2 _ > r A;-Riti dA;Riti (a Q\ 

\ U A,A/a;mi M I U k,MA^2,V2l - a nuvi H )a,R 2 ,T2^^,U2 — A W^i,fi a H2,V2 ■ l D - y J 

These states will be orthogonal as long as the eigenvalues differ A Ml 7^ A M2 . In the 
half-BPS case, V — I and VQV ~ (Dimi?) -1 while the states are labeled by R and so all 
eigenvalues are distinct. Also for all cases A = [n, m] with m > 2 that we checked explicitly 
we did not find any degeneracy among the eigenvalues. The only exception seems to be 
the sector A = [n, 1]. With only one X operator in the traces symmetrization doesn't 
act and VQV ~ (Dimi?) -1 again. But now there is a possible multiplicity in r and two 
distinct states with the same R will be eigenstates with the same eigenvalue (Dimi?) -1 . 
This happens first at A = [5, 1] where we find operators Oh=[3,2,i],t=1j @r={3,2,i],t=2- If h 
is true that this pattern is general, namely that the cases of degenerate eigenvalues for 
fixed A only arise with A of the form [n, 1], then for most A the eigenvalues of VQV can 
provide a way to label the orthogonal eigenvectors of the 2-point function. 

Even though we can write the orthogonal basis formally as (16.81) . finding the eigenvec- 
tors B^' Rt in practice is of course a very difficult task. We provide an example result in 
the Appendix IC. 1.31 for the simplest nontrivial case A = [2,2]. In general they come out 
to be complicated expressions involving N. We will not have much more to say about the 
precise form of the eigenvectors. 
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6.4 Towards a complete set of labels for the BPS operators 

We will now describe a matrix related to VQV whose diagonalisation should be simpler 
and which may provide a more practical way of arriving at a labelling of the eigenstates 
of VQV . Consider the large N expansion of Q : 

S = l+^Gi + o(^). (6.10) 



N \N 2 



{yx)R 2 - -oruR* — - A (o- 1 - 1 -; 



with S[2] being the sum over all permutations with cycle structure [2]. Note that XRl ^ l2] 
is just a number and not an iV-dependent function. The leading terms of VQV are then 

VQV = V + ^VGiV + O (J^j . (6.12) 

V and VQ\V commute due to V 2 = V and they can be simultaneously diagonalized. 
This implies that the eigenstates of VQV in the limit iV — > oo will approach those of 
simultaneous eigenstates of V and VQ\V . The BPS states will correspond to the subspace 
of V eigenvalue 1 and their eigenvalues will be 

A , = 1 + ^ + °(^) ( 6 - 13 ) 

where the A M are the eigenvalues of VQ{V. 

It is conceivable that, for fixed A, the eigenvalues of VQV do not cross when plotted 
as a function of N. Indeed, this is the case in the examples we have studied concretely. 
In that case, the eigenvalues of VQ\V would be adequate as a way to label the states 
and identify them systematically with giant gravitons and their excitations in the regime 
n ~ N, and with LLM-like geometries in the regime n ~ N 2 . This deserves further 
investigation. 

Progress on the dictionary between space-time giant graviton and LLM-type states 
will also require the computation of appropriately normalized three-point functions. For 
example we may expect the three point function of an operator corresponding to a giant 
graviton along with some spacetime excitations to be of order one when the other operators 
in the 3-point function are the giant in question and the small operator for the space-time 
excitation. The techniques of [IB] using correlators in connection with the asymptotics of 
SUGRA solutions should be useful in relating correlators of eigenstates of VQV to bulk 
spacetime. 
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7 Summary and Outlook 



We have described an elegant method for a refined counting of quarter and bosonic eighth- 
BPS states, where states are organized according to representations of U(2) or U(3). The 
method extends to general U(M). It relies on the construction of an element P living in 
the group algebra of S n . Exploiting the Schur-Weyl duality relation between represen- 
tations of U (M) and symmetric groups, the characters of P in different representations 
of S n provide the multiplicities of different U(M) representations A among BPS states 
constructed from n matrices, chosen from X 1; X 2 , ■ ■ -Xm- This result is given in ( 13 .31)) . 
The coefficients of different conjugacy classes in the expansion of P match in simple cases 
known integer sequences, and generalize them (see Appendix |E]). Recent results on the 
generating functions for BPS states [131 ESj in terms of simple harmonic oscillator states 
in M dimensions, after a judicious change of variables (13.421) . provide the generating 
functions for these coefficients of conjugacy classes in P. 

The fine-structure of P encoded in a map p : C(S n ) —> C(S n ) with matrix elements 
Pp, a lead to the construction of a linear operator V acting on the Hilbert space of the free 
M = 4 SYM theory (i.e the theory at g\ M = 0) which projects to the symmetrized traces 
which are BPS states in the planar limit, annihilated by the one-loop dilatation operator 
[32] . The matrix elements of V are given in terms of a representation theoretic basis 
for the free theory, which diagonalizes the free two-point functions at finite N [161 02] • 
The free basis is a Fourier transform of the trace basis. The expression for V in this 
Fourier basis is given in terms of Clebsch-Gordan coefficients and matrix elements of the 
symmetric group ( 14. lip . 

BPS states to all orders in the large N expansion are given in (I4.14p as the image of 
QV, where Q has simple matrix elements in the Fourier basis, given in terms of inverse 
dimensions of representations of U(N). The two-point functions of the BPS states are 
given by the matrix VQV (14.261) . The matrix V also has the virtue that its kernel is 
spanned by the descendants, which are not annihilated by 7^2- As a result it allows a 
manifestly finite A" construction of BPS states (15. 1 8 p and their two-point functions (15. 19[) 
ViQVi. The eigenstates of VjQVi can be used to construct higher point correlators in a 
manifestly finite A" setting. Analyticity makes it possible to reconstruct eigenvalues and, 
conjecturally, eigenvectors of the ViQVi from the matrix VQV . Non-renormalization 
theorems allow us to argue that the eigenvalues of the matrix VQV contain information 
about gauge theory operators for giant gravitons [11] and their spacetime excitations and 
LLM-like geometries [HI H5l 061 07]. The use of the stringy exclusion principle allows 
some progress in this direction. An important open problem is to find a labelling of the 
eigenvectors of VQV, which will allow the identification of different types of semiclassical 
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objects : giant gravitons and their spacetime excitations, LLM geometries and their 
fluctuations. We discussed some avenues in this direction in section |6j The importance 
of the matrix VQV cannot be overstated^- 

The route leading from P, which furnishes results in counting, to the matrix V and the 
two-point functions VQV is very interesting. It seems to capture a sort of categorification 
[4"9] . where numbers (multiplicities) are promoted to explicit constructions of quantum 
states in string theory, in a manner that can make semiclassical limits and space-time 
physics transparent. It will be interesting to find out if a similar route can be followed 
in more general systems in string theory, to lead from results in counting to explicit 
operators, states and semi-classical limits. 

We end by mentioning some additional avenues for further research. 

• Generalization of P, V and V1QV1 to the case of the most general eighth-BPS states 
with global symmetry group £7(3 1 2). Other generalizations to consider are for the 
SU(N),0(N),Sp(N) gauge groups. 

• Developing the identification of eighth-BPS giant graviton states further and finding 
operators corresponding to strings attached to giants, generalizing the use of the 
Young diagram description to characterize strings attached to half-BPS giants [50] 

EH E21 E31 E3 

• The matrix VQV can be described in other bases which diagonalize the free 2-point 
functions. The transformation between the restricted Schur basis [SHUSHES] and the 
Fourier basis is described in [5?]. A similar transformation can be used to connect 
to the Brauer algebra basis [191 EH]- This should shed light on the construction of 
[59] which leads to an elegant description of a subset of quarter-BPS states in terms 
of projectors in Brauer algebras. This could provide another angle on the problem 
of finding a neat labelling of the eigenstates of VQV which connects to spacetime 
physics. 

• By exploiting a result of [60] on the one-loop dilatation operator for the U (2) sector, 
the structure of Clebsch-Gordan coefficients and symmetric group elements which 
we used to describe V can also be used to describe the action of the one-loop 
dilatation operator on the finite N Fourier basis in a manifestly U (2) covariant way 
(see Appendix |B]). The study of this mixing in the Fourier basis was initiated in 
[61j . The equation (IB.lOj) should shed light on the one- loop mixing problem of the 
descendants, which are in Ker(P). The mixing problem at planar level and beyond 

4 We will in fact confess that an earlier title of this paper which did not make the final cut was : "The 
Matrix : VQV? 
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has been a very active subject of research, notably in connection with integrability. 
A few references, by no means complete, giving the flavour of the subject are [26| 

E21 m EH EH E3 ESI m\ 

• A better understanding of the quarter and eighth-BPS sector is expected to have 
implications for black hole physics EH]- The identification of gauge theory 
operators, involving impurities added to the eigenstates of VQV, corresponding 
to black hole solutions will be interesting to explore. In [70], there was evidence for 
commuting variables in the counting of sixteenth BPS states, with analogies to the 
eighth BPS sector. This suggests something akin to the V operator could be useful 
for the sixteenth BPS case. 
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A Further properties of V and related operators 

The operators V represents the action of the symmetrization operation in the basis 
@a,m a ,r,t which diagonalizes the free two-point function. We describe the properties 
of V and of some closely related projectors which are used in counting of BPS operators. 

A.l Glossary of operators related to V 

¥ : Introduced in (I3.30p . A universal element in C(S n ). Key counting result (13.311) . 

P p : Introduced in (13.261) . It is a projector in C(S n ) which depends on the partition 
p. Relation to IP given in (13.301) . 

¥ a : Introduced in (14.11) . also a projecctor in C(S n ). It is a slightly refined version of 
P p with a specified embedding of p. 
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p : Introduced in (I4.2p . It is a linear map from C(S n ) to C(S n ) with pp j0l as the matrix 
elements of the map. p^^'] is the Fourier transform of pp t0l . 

PW : Introduced in (IXlgjh The finite A" version of the universal element P captures in 
symmetric group language the fact that the chiral ring counting is equivalent to counting 
states of N particles in a simple harmonic oscillator. 

V : Introduced in (14.111) . It is a matrix built from p, using Clebsch-Gordan coeffi- 
cients of S n . Key property: it is a projector implementing trace-symmetrization in the 
Fourier-basis. 



V\ a ] : Introduced in (lA.lOj) . A summand that appears in the expression for V. Also 



a 



projector in the Fourier basis which symmetrizes only the given trace structure [a]. 

A.2 Poperties of p : C(S n ) -> C(S n ) 
The equation ¥ 2 a = F a implies that 

^ Pp,aP^T,ot = Pra (A.T 

P 



We can convert the equation flA.lj) into an equation in terms of the Fourier transformed 
quantity. 

RlPiqi _ \ d,S 1 ds 2 Si h jl S 2 12 32 qS lt S2\Rl,T oS 1 ,S 2 ]Rl,T /A n\ 

"-R2PK72 / j "-R2 P2 m Pr 2 m q 2 ° jl 32 91 »l *2 PI \ " / 

Sl,S 2 ,tl,l2 ,31 ,32 

Figure [3] gives a diagrammatic expression of the equation. 
The coefficients pp, a have the symmetry 

PP,a — P7 / 37- 1 ,7a7- 1 (A. 3) 

Hence ^a-sP/V /3 ® a defines an element of C(S n ) <g) C(S n ) which commutes with the 
diagonal C(S n ), and contains all the information about the counting Xa{^o.) of irreps A 
among the symmetrized traces of trace structure detrmined by the cycles of a. Such 
elements were considered in [35] in order to resolve the r multiplicity of symmetric group 
Clebsch-Gordans. 

Pp~\a = P/3,a (A.4) 
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Using -Dg(/5 -1 ) = Dfi(P). This means that 

Psl J i = PsH (A.5) 

We also have 

P/8,a = P 7 /3,a = P/3 7 ,a for 7 G G(o) (A.6) 

where is defined in section |4~T1 
A. 3 Projectors V\ a ] 

The matrix V in ( 14. lip can alternatively be written as: 

(P)ii£S = £ A*(PJ (ff,U«)St ") (D% h (a)S% * *• ") . (A.7) 

This is seen by writing out = J2 a D^ Am iJJ? a )Dl ik2 (a) and then performing the 

sum over the representation S. Figure 0] expresses the RHS of (1A.7|) in a diagramatic 
form. 

The first thing to note is that the argument of the a sum in (I A. 7ft is invariant under 
conjugation a — > / ya , ~f~ 1 . That is because 

SgOPWO = A^TPaT" 1 ) = Av(Pa)Ai(7)^y( 7 ) (A.8) 

and 

(A.9) 
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and so the resulting -D A (7) matrices multiply to the identity. This invariance means that 
the sum over a can be really seen as a sum over the conjugacy classes [a] each coming 
with a factor of n!/|Sym(a)|. It is then convenient to define a new matrix V[ a y. 

where the right hand side can be evaluated with any permutation a in the same conjugacy 
class, and thus V[ a ] only depends on [a]. The full matrix V is then just a sum over 
conjugacy classes or partitions: 



[a]\—n 



The newly introduced matrix V[ a ] can be seen as P Q or P p projector in C(S n ) appro- 
priately transformed to the |A, R, r) basis in the same way that V is a transformed P. 
The most remarkable fact is that V[ a ] itself is a projector with each [a] projecting to a 
distinct subspace! That is 

V[ a \P\0\='P[ a ]8[ a ^\p\. (A. 12) 

This can be shown by evaluating 

J2d R (D^ h (a)S\\ A l T ) (d^S**^ t ) = Df^aa^a^) (A.13) 
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and then using this result to take the product: 



E,„ \A,Ri,ti /-r> \A,R,r 



s A,R,t 

!,T2 

R.T 



R 2 R2 A, t 2 
k 2 h j 



^MWiZ. | Sym(a) | 



(A.14) 

In the third line we used the fact that with the delta function present we could replace 
(3 = a- x aa and that = P Q . 

It follows from (jXl2l) and flATTI) that V itself is a projector 

V 2 = V (A.15) 

which is a nice consistency check. 

It is interesting to ask what is the significance of the individual projectors V[ a ], that 
is, what operation on the traces do they correspond to. More precisely, let us define a 
linear operation on the Hilbert space 

symm w [O AMa , r , t ] = £ (T^ffi^O^^ (A.16) 

and ask what 

symm [CT] [tr n (X s a)] (A. 17) 

evaluates to. From the derivation (14. 9 j) and (I4.10p it can be seen that in order to get P^j 
on the RHS rather than "P the sum over a has to run only over the conjugacy class [a]. 
That implies that the following action is needed on the trace basis instead of ( 14. 31) : 

symm [CT] [tr n (X s a)] = ^ P/3,atr n (X /3(( j ) a). (A.18) 

The interpretation is clear: symm^ and its representation on the Fourier basis V[ a ] acts 
by symmetrizing only the multitraces with the given trace structure [a]. All traces with 
different structures are annihilated. This explains (IA.12I) : if [a] and [/3] are different, 
there are no traces which survive. Also (lA.llj) has a straightforward interpretation that 
full symmetrization is performed by symmetrizing all trace structures. 
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B One-loop dilatation operator in the Fourier basis 



In this section we derive an expression for the one-loop dilatation operator % 2 i n the 
quarter-BPS sector acting on the Fourier basis. We use the same methods that led us 
to the V matrix for symmetrization. The goal is to find the matrix (H2) \ 'mv" 1 '/? 1 ' 1 " 1 sucri 



>A2,M' An ,R2,T2 

that 

^Ai,M Al ,Bi,-n 

" / A 1 ,M Ai ,R 1 ,t 1 
Ai,M Al ,R u n 

and to express it in terms of symmetric group objects. The matrix (%) could allow us to 
go beyond the analysis of BPS states (which correspond to its zero-eigenvalue eigenvectors) 
and to calculate the anomalous scaling dimensions of the non-BPS operators. 

Inspecting the derivation of V in (14. 9p . fl4.10p we find that the main ingredient there 
is the action of the linear operator of interest on the basis. In the case of V that was 
the symmetrization operator 

symm[tr n (X s a)] = ^ p^tr^X^q). (B.2) 

The V matrix is then a transformation of this action to the |A, Ma, R, t) basis. 

The one-loop dilatation operator in the U(3) sector acting on multitrace operators can 
be written as 

U 2 = --tr([X l ,X J ][X u X J ]) (B.3) 

with 

(X a )i ee 9 . . (B.4) 
This can be translated to the Og a basis as [60] 



% 2 tr n (X 3 a) = y^tr ra (X Pi .( S ) (ia(j))'a ). (B.5) 

Here is an element of the symmetric group algebra 

Pij = l-(ij), (B.6) 

with the action on X^ defined as Xp t .p) = X^ — X^j)^). The term (ia(j))' is a permutation 
interchanging i and a(j) which multiplies a, except for when i = a(j). In that case it 
should be understood as a factor of N: 

r . v _ /(#) if i±h m v x 

M =< M . (B.7) 
I A* if i — J. 
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It is an element of C(SVi) when N is viewed as N £ C. We can view flB.5j) as expressing 
the fact that H2 is determined by an element 

«®(i«0'))'®^i (B-8) 

in C(S' n )(8)C(5' n )®C(S' n ), which is an analog of the operator p in C(S n )(&C(S n ) introduced 
before. 

Now that we have expressed the I-L2 action in terms of C(S n ) quantities, we can 
transform it to the Fourier basis. Repeating the derivation along the lines of 04. 9p . (14.10]) : 

-H 2 O a , Ma , r , t = ^EE^ ?m T ^(«)C , f,M A , m tr n (X Pfei(5) (ka(l))'a ) 

k^l a, a 



k^l a,a 



E y/^C^- mi D^((ka(l)ya)S * f^A,^,^ 

Ai,M Al ,Ki,ri 



E, y/diidih 
U *,Ma,Ri,ti -j 



2)A,M a ,R,t ^Ai,M Al ,-Ri,7i- 
Ai,M Ai ,_Ri,ti 



We have thus calculated the transformation matrix to be: 



(B.9) 



Ai,M Al ,fli,ri _ a/ d^dfy 



x E E ^U(^) (^((MO)'") s t J ir) ^ 1 

k^l aS5„ 



R 2 Ai,r 2 
42 J2 T12 



(B.10) 



This expression should be compared to (IA.7I) . One can see the D A (P k i) appearing in the 
same place where we had _D A (P a ) before, because that is the permutation acting on a in 
the trace basis. A new ingredient here is the (ka(l))' operator acting on a, which was 
absent in the symmetrization action. 
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C Examples of BPS operators 



C.l A =[2,2] 

In this section we will provide an example of explicit calculation of BPS operators in 
the representation A = [2,2]. Since the operators throughout the section will all have 
the same A, also M\ will be the highest weight state and r = 1 because there is no 
multiplicity, we will abbreviate 

Or = O 'A=[2,2],A/ A =HWS,ii,T=l- (C-l) 



C.l.l Large N operators 



First we will assume N > 3, but will keep all 1/N corrections. Then in the next section 
we will also explain what happens for lower N. 

We start with the basis in the free theory. Using (12.41) the operators are: 



Cr=[3,i] 



1 



4 v / 3A r2 
1 

2VQN 2 
1 



(-tr(X i X i X i X i ) + tr(X i X i )tr(X i X- J ') + tr^X^tr^trpf-*)) 
(ti'(.Y,.V'.V ; .V r ) + trtXiX^trpntrpS^')) 



4V3N 



2 (-tr(X i X i X j X-'') - ^(XiX^tviX^^ + tr(X i X j )tr(A: i )tr(A: J ')) . 

(C.2) 



Here for compactness we are using upper-lower index notation meaning XiX % = X\Xi — 
X 2 Xi, which is a singlet in the quarter-BPS sector. 

The symmetrization matrix V can be evaluated using (14. lip or (IA.7[) . In the A = [2, 2] 

sector we get: 



V 



( 3 jl _i \ 

4 2^ 4 » 

1 1 _J_ 

\ 4 2V2 4 / 



(C.3) 



It can be verified that transforming the free operator basis with this matrix indeed sym- 
metrizes terms in each trace. In this example that just means that the tr(XiXjX l X^) is 
dropped. 



( °[3,1] °[2,2] ^2,1,1] ) = ( ^[3,1] °%2] C[2,l,l] ) 



/ 3 _±_ _1 \ 

'4 2^ 4 » 

1 1 _J_ 

2^ 2 2^ 

_1 _J_ 3 

4 2^ 4 



(C.4) 
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o. 



[3,1] = (tr(X^)tr(X^) +tr(X,X,)tr(X*)tr(^)) 



Of 2>2] = (tr(X^.)tr(XOtr(X^)) 
1 

• 1] ~ aVsn 2 



(C.5) 



(-tr^X^trpTX^) +tr(X i X i )tr(X i )tr(X j )) . 



Clearly, there are only 2 independent symmetric operators. This corresponds to the fact 
that at large N there are only 2 BPS operators when the coupling is turned on, compared 
to 3 operators in the free theory. It is also worth noting that the null eigenvector 



-i 



(C.6) 



corresponding to |C[3,i] — ^C[2,2] + |C[2,i 2 ] is the descendant operator. This is a general 
fact that the descendants can be characterized as Ker(P). 

We can write down the precise expressions for the BPS operators using (I4.14p : 



(C.7) 



Run 



with the diagonal matrix 
Q = diag 



N n d R \ 
n\ Dim/? J 



TV 3 



N 2 - 1 



\ 



\ 





1 

N-2 



(C. 



/ 



Let us arbitrarily pick OWS and Oj^^ as the two linearly independent operators. In 
terms of the free operators we get 



U [3,l] 



/-OBPS 



AT3 



3 O 



[3,1] 



1 0[ 2 ,2] 1 C[2,l,l] 



iV 2 -lV4iV + 2 2^2 N AN -2 



N 3 



N 2 



1 0[3,1] 1 0[ 2 ,2] 3 C[2,l,l] 

'In + 2 2^/2 n In -2 



(C.9) 
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Explicitly in terms of products of traces they are 



O 



,BPS 




[3,1] 




+N(N - l)tr(X i X i )tr(X i X i ) + (N 2 - 2N - 2)tr(X i X j )tr(X 4 )tr(X J ')) , 



o 



,BPS 



[2,1,1] 




-N(N + l)tr(JC i X j )tr(X i JC i ) + (iV 2 + 2X - 2)tr(X i X j )tr(X i )tr(X i )) . 

(CIO) 



To the leading order in TV the operators Oj, are equal to the symmetric combinations 
0%, but there are further 1/N corrections which make the state exactly annihilated by the 
one-loop dilatation operator. Note how the antisymmetric combinations tr(XjXjX J X- J ) 
appear in the exact operator, although only in the subleading-X order. 

C.1.2 Finite N cutoff 

The expressions for BPS operators (IC9p . fIC.lOp are valid whenever N > 3. Clearly 
something goes wrong when N = 2 as indicated by N — 2 in the denominator - this 
is where the finite- X cutoff takes effect. We will demonstrate here how it is dealt with 
according to the general procedure described in Section It turns out that at N = 2 
there is only one BPS operator remaining. 

The cutoff is best seen in the free operator basis: Oa,a/ a ,r,t becomes whenever Young 
diagram R has more than N rows. In our case one of the three free operators drops out: 



In the expressions for BPS operators exactly these free operators which vanish are accom- 
panied by in denominator because of Dimi? = 0. 

It is not immediately obvious how then to find correct BPS combinations. One could 
try to just naively cut off the terms in (10.9H and write 



C[2,i,i] = 0. 



(C.ll) 



Note that implies a relationship between products of traces, using (1C.2|) : 



tr(XiX l XjX 3 ) - tr(XiXj)tr(X l X 3 ) + tr(X i X i )tr(X*)tr(X') = 0. 



(C.12) 
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however, this is not correct. If it were true, it would seem that we are still left with two 
BPS operators, however, it turns out that even at N = 2 neither of them is annihilated 
by one-loop dilatation operator! 

What we have to do instead is to pick linear combinations of C^ PS which do not depend 
on these Or that are set to zero. These linear combinations will be "protected" from what 
happens to the terms with Dimi? = in the denominator. This procedure corresponds 
exactly to picking vectors in the intersection space lm(Vi) described in Section [5j where it 
is shown to give all the BPS operators. In our example we have to take the combination: 



qBPS;(N=2) 



U [3,l] 



1 /^BPS 

3 U [2,i,i] 




(C.13) 



This is the single remaining BPS operator at iV = 2. To be more precise, in the language 
of Section [51 we find the projector Vj for the intersection of Im(Z(jv)) and lm(V) to be 



V\ 



(JV=2) 



/ 


2 
3 


V2 
3 
1 





\ 




V2 









3 


3 















) 



(C.14) 



Then we get the BPS operators by acting with {GVi) on the full Hilbert space, which 
gives precisely the one-dimensional space spanned by the operator ( 1C.13I) . 

There is a related subtlety worth pointing out. It is true that counting BPS operators 
is the same as counting states in the chiral ring, where the matrices Xi can be taken to 
be diagonal. One may wonder whether this is the same as counting independent products 
of traces, where the terms in each trace are symmetrized, that is, the 0\. Rr operators. 
It turns out that this is not the case when the finite-iV cutoff applies. In our present 
example at N = 2 there are still two independent symmetric operators: tr(XjXj)tr(X i X- 7 ) 
and tr(XjX J )tr(X*)tr(X- ? ). However, one has to eliminate one linear combination - the 
one that is equal to tx{XiX l XjX^) by the relationship ( IC.12j) . This leaves us with one 
operator. 



C.1.3 Orthogonal BPS basis and VQV eigenvalues 

In order to find an orthogonal basis for the BPS operators we have to look for the eigen- 
vectors of the two-point function matrix on the BPS operator space VQV . In the case of 
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our present example A = [2,2] we can calculate the matrix to be: 

/ 



N 

VQV 



(N 2 - 1)(N 2 -4) 



V 



3N 2 -4N-2 


N 2 -2N-2 


N 2 +2 


4 


2V2 


4 


N 2 -2N-2 


N 2 -2 


7V 2 -27V-2 


2V2 


2 


2V2 


N 2 +2 


N 2 ~2N~2 


3N 2 -4N-2 


4 


2V2 


4 



(C.15) 



Its characteristic polynomial is 



det(VgV - XI) = (Ar2 _ 4) ^ 2 _ 1)2 {-(N 2 - 4)(iV 2 - 1) 2 A 2 + 2iV 2 (iV 2 - 1) 2 A - iV 6 ) 

(C.16) 

and we find eigenvalues 

^ ( ^-4)(^-i) (jV-iTVSvr+i). (c i?) 

A 3 = 0. 

The eigenvectors with zero eigenvalues are the same as those of V and span the non-BPS 
combinations. We are interested in the other two eigenvectors 

BPS -3N T 2V2WTT B p S -V2(N+l±V2WTl) BPS N- 2 BPS 

U (i,2) - ]y U [3,l] + jf U [2,2] + —^—^[2,1,1] 



A^ 2 



N 2 - 1 



-3AT T 2V2N 2 + 1)^^ - V2(N + 1 ± V2iV 2 + 1)^M + [211] j 

(C.18) 



These two operators provide us with an orthogonal basis for the BPS sector. The normal- 
ization that we picked so far is arbitrary, except that we adjusted it to allow for a nice 
limit to A" = 2 as we will see shortly. 

C.1.4 N = 2 limit revisited 

As discussed already, at A^ = 2 there is only one BPS operator in the A = [2, 2] rep- 
resentation given by the combination flC.131) . On the other hand, the operators 0^ s , 
0^ s provide an A^-dependent basis of orthogonal states from a calculation done assum- 
ing N > n. According to the discussion in Section 16.11 we can in fact take the limit of 
the orthogonal basis to A" — > 2 and find the one remaining BPS state as the surviving 
eigenvector with finite eigenvalue. We will see an example of that here. 

Let us simply take the expressions for eigenvalues and eigenvectors flC.17j) . (1C.18I) and 
calculate the limit A" — > 2. We find 

\{ N=2) = -, Af =2) ~-^-^oo (C.19) 
1 9 2 N-2 v ; 
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and 

( '\\) ~ —*^[3,lj - -*v ^[2,2j -r - 



BPS;(7V- 2 ) = _ A ^/2 0[2 2] + 4^ 

4 3 (C20) 

~ 3' 



U (2) ~ o U [2,l,l]- 



The operator 0[2,i,i] is identically zero in this case and so can be dropped from the 
operators. The first eigenvector 0?^ S 2 \ which has a finite eigenvalue, is indeed the 
right BPS state (1C.13I) up to an overall normalization. The second eigenvector, with 
divergent eigenvalue, disappears from the Hilbert space. This reflects a pattern that we 
expect to be general. 

C.1.5 N — > oo limit revisited 

Another interesting limit of the orthogonal BPS operators that we can take is iV — > oo. 
We know already that at large N the leading term of the BPS operators is always a sym- 
metrized trace. So a priori we can take any orthogonal combination of tr(X i Xj)tr(X l X : >) 
and tr(XiXj)tr(X' l )tT(X^) to be the large-iV basis. What is interesting about taking the 
limit of ( 1C.18I) is that it provides us with a unique preferred basis. We get 

bpS;(*=oo) = ( _ 3 T 2 ^2)O m + (-V2 =p 2)0 [2l 2] + 0[2,1,1] 

2±a/2 1±a/2 ( C - 21 ) 

7 -tTiXiXjMxw) - -^tTiXiXjMX'Mx*). 




2V3N 2 1 3 

Both eigenvalues approach 1, because VQV — > V. These states can also be found as the 
eigenstates of VGiP, as discussed in Section 16.41 The Q\ matrix is in this case: 



(C.22) 



C.2 A =[3, 2] 

Here we work out the calculations for the A = [3, 2] sector. In this section we again 
abbreviate the operators: 

Or,t = C>A=[3,2],M a =HWS,K,t ^ 
Or = C , A=[3,2],M A =HWS,iJ,r=l 

with the U (3) label Ma in the highest weight state. 
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The free operators: 



o R=m = 2V ^ N5/2 (-tr(X 1 X i X J X,X^) - ^(X^X^MX,) 

+tr(X 1 X i X j )tr{X i X j ) + -tr(X 1 X l X i )tr(X i )tr(X J ) 

3 

+hi(X i X j )tv{X i Xi)tr(X l ) + UriXXjMX^triXiMXr) 
° R=m = 4^2X5/2 (-to&iXiX***') ~ |tr(X*K%tf)*(*i) 



--tr^iJCiXjOtr^'Jtr^') 



2 

-itr^XjOtr^X^tr^i) - ^tr^X^tr^^tr^^tr^)^ 

R=[3>l2] , r=1 = 4 ^ v5/2 (tr^XX'X^) - 2tr(X 1 XX J )tr(X 4 )tr(X^') 
+tr(X i X j )tr(X l X J )tr(X 1 )) 



0JM*M] = A J N5/2 (-1x(X 1 X i X i X j XS) + ^tr(XX i X,X^')tr(X 1 ) 

-ItTiXiXiXjMx'Mx*) 

-^tr(X i X J )tr(X i X^)tr(X 1 ) + ^tr(XX i )tr(X 4 )tr(X^)tr(X 1 ) 

° R=[2 ' 13] = UiXxXipXjX') - hnixjPXjXtMXi) 

+ti(X l X i X j )ti(X i X j ) - -tr(XiX i X i )tr(X i )tr(X i ) 
-^tr^XjOtr^X^tr^) + ^tr(XX i )tr(X t )tr(X^)tr(X 1 )^) 
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The symmetrization matrix: 
/ 



V 



2 
3 

6 

1 
2 

5^ 
1_ 

6V$ 

J_ 

15 



6 

_r_ 

12 

1 



1 

hs/2 
1 



2V10 
1 

12 

1 

6^ 



2V10 
_9_ 
10 


i 



2V10 
1_ 

5V? 



2 


1 


1 


\ 


5^ 




15 


1 


i 


1 




VIE 


12 


6^ 







1 


_ 1 




2^ 






3 


1 


2 




5 


VTB~ 


5^ 




1 


7 


V5 




Vv> 


12 


6 




2 




2 


) 


5^ 


6 


3 



(C.25) 



which has rank 4, corresponding to four BPS operators in this sector. The Q matrix: 



N 4 



(N 2 -1)(N 2 -A) 



N-2 
7V+3 









The BPS operators are again given by 



O 



BPS 

R,T 



E 





N-2 
N 






Rl,Tl t 







1 
1 










7V+2 
TV 











7V+2 
7V-3 



\ 



(C.26) 



(C.27) 



VQV 



The resulting combinations are best seen by inspecting matrices V and Q and not worth 
writing out explicitly. 

The two-point function matrix on the BPS states will be 
N 3 



(N 2 -1)(N 2 

/ KW_ 3 -37A r2 +llA r +36 



15 



-4) (TV 2 -9) 

57V 3 - 237V 2 +47V+54 

77V 3 -167V 2 -377V+72 
12 



7V 3 -37V 2 -6 


27V 3 -47V 2 +47V-18 


7V 3 -37V 2 -47V- 


-IS 


7V 3 +117V 




5^ 


6V$ 




15 


7V 3 -157V+12 


7V 3 +7V 2 -137V+9 


7V 3 -197V 




7V 3 +37V 2 -47V+18 


2VW 


VTE 


12 




6^ 


97V 3 - 757V 


67V 2 -18 


7V 3 -157V-12 




7V 3 +37V 2 +6 


10 




2^ 




5V2 




37V 3 -197V 


N 3 -N 2 -13N 


-9 


27V 3 +47V 2 +47V+18 




5 


VIE 




5^ 






77V 3 +167V 2 -377V- 


-72 


57V 3 +237V 2 +47V-54 






12 




6^ 










107V 3 +377V 2 +117V-36 










15 



The characteristic polynomial of this matrix is: 

\2 

det(VgV - AI) 



(N 2 -9)(N 2 -4) 3 (N 2 - iy 

\3 \3 



((iV 2 -9)(iV 2 



— 4iV 4 (iV 2 - 4) 3 (iV 2 - l) 3 A 3 + 3iV e (2iV 4 - 3iV 2 + 4)(vV 2 
-2N 10 (2N A - 3N 2 + 4)(iV 2 - 1) A + iV 16 ) 



(C.28) 

-4) 3 (iV 2 - 1) 4 A 4 
4)(iV 2 - 1) 2 A 2 



(C.29) 
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The A 2 prefactor indicates two zero-eigenvalue states corresponding to the kernel of 
V. The remaining fourth-degree characteristic polynomial equation can not be solved 
explicitly, but it is useful for inspecting what happens at finite values of N. When N 
approaches critical values, some eigenvalues of VQV run off to infinity while the number 
of eigenvalues that remain finite is equal to the number of surviving BPS operators. One 
can see that (1C.29j) reduces to a third degree polynomial when N — > 3 and to a first degree 
polynomial when N — > 2. We can conclude the following multiplicities of BPS operators 
depending on N: 



N 


# of BPS states 


> 3 


4 


3 


3 


2 


1 


1 






Finally, we find the orthogonal BPS states in the limit of iV — > oo. That is, the 
preferred basis of symmetrized traces. We calculate the 0(1/N) term of Q: 



I 



N' 



( -5 





V o 










+ 0(l/iV 2 ), 

\ 







2 

5 / 



(C.30) 



(C.31) 



Then the eigenvectors of 



VQ{P = V + jjVQ{P + 0(1/N 2 



(C.32) 



approach the eigenvectors of VQ{P as N — > oo. Since V and VQ{P commute, they can 
be simultaneously diagonalized. That means that we can resolve the degeneracy of the 
symmetrized traces which all have V eigenvalue +1, by finding combinations of them 
which are at the same time eigenstates of VQ\V, but with different eigenvalues. We find 
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these eigenstates reexpressed in terms of multitraces to be: 

o bpS;(j\t=«0 = iotr(XiX i X J )tr(X^ J ) + (-Vl9 + 3v / lT)tr(X 1 X i X j )tr(X 4 )tr(X J ) 

+ ^(3^ + v / lT)tr(X i X i )tr(X' i X J )tr(X 1 ) + -(-3 + v / 209)tr(X i X j )tr(X i )tr(X J ')tr(X 1 ), 

hps ; (jv=oo) = 10 tr(X 1 X i X i )tr(X^ J ) + (-v 7 !^ - 3v / TT)tr(X 1 X i X J )tr(X i )tr(X J ) 

+ ^(3^ - v / Tl)tr(X i X.,)tr(X i X J ')tr(X 1 ) + -(-3 - v^MXiX^tr^tr^trpfj), 

hps ; (jv=oo) = xotr^iXiX^tr^^') + (>/l9 + 3v / TT)tr(X 1 X 4 X j )tr(X l )tr(X J ) 

+ -(-3^19 + v / IT)tr(X i X i )tr(X i X''')tr(Xi) + -(-3 - v / 209)tr(X i X i )tr(X i )tr(X i )tr(X 1 ), 

bps;(n=oo) = 10tr (X 1 X i X j )tr(X i X J ) + (>/l9 - 3v / lT)tr(X 1 X i X J )tr(X 4 )tr(X J ) 

+ -(-3\/l9 - v / TT)tr(X i X j )tr(X i X J ')tr(X 1 ) + -(-3 + v / 209)tr(X i X J )tr(X i )tr(X i )tr(X 1 ). 
2 2 

(C.33) 

C.3 Computational remarks 

The main challenge in the calculations was getting the V matrix. In our implementa- 
tion we used Mathematica to evaluate the expression ( IA.7j) . It involves three ingredients: 



calculation of P a , the representation matrices D^(a) and the Clebsch-Gordan (CG) coef- 
ficients ^ fc A ' T . We evaluated the algebra element P a directly using the definition ( 14. ip . 
The representation matrices were calculated using the method described in Section 7-7 of 

Finally, the CG coefficients S i • k ' T were calculated using the recursion relationships 
according to the method described in Section 7-14 of [ID]. This is the most computa- 
tionally intensive procedure. We also found that the relevant equations (7-226), (7-233) 
needed a small clarification for dealing with multiplicities, which is worth pointing out 
here. The coefficients Kj^f used in the recursion relationship have, in fact, two extra 
multiplicity indexes. The equation (7-226) should read as: 

s s z 5 r = E K " r V) sj T \ {cm) 

where r' runs over the multiplicity associated with (7$, a a , j3 e ). Then Kj " /' T (r) is the 
r'th solution of the equation (7-233) which should read as: 

t>-7 ol P, t' o7i «c Pg, t' i ra fP _ nf\ _i_ t>-7 a P, r' 07, a c B h , r' ra 8 
Iy i c g °jt dr hs \J cdJ gh Jij) tn ic/i °jt dr gs J cdUhg 

I jsl a P, r' c7i <*d Pg, t' at fP , rsl a P, r' Q7i "d Ph, r' a P _ jsl a B, t' qjj a c Pg, r' 7 

xy i d g °jt cr hs UdcJ gh "i ^ i d h °jt cr gs ildc'dhg ~ xy j c g °it dr hs ^ij^ 

(C.35) 
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with sums over i, j, c, d, g, h, t, r, s and the extra multiplicity r'. 

The Mathematica code is available from the authors upon request. 



D Multiplicity tables 

In this appendix we give some counting information. 

First we write out the evaluated universal element as defined in (I3.30p . We organize 
the terms in increasing n. 

F(S 2 ) = I (2E [2] + 2Z m ) , 

P(5 , 3) = i(2E [3] +3E [2)1] + 5E [1)1)1] ), 

PC&i) = (3S[4j + 3S[ 3)1 ] + 7E[ 2) 2] + TEp.i.i] + 15£[i, 1,1,1]) , 

P(S 5 ) = ^ (2E [5] + 4S [4 ,i] + 5E [3)2] + 7E [3 ,i,i] + 12E [2 , 2>1] + 20E [2)ltl>1] + 52E [1 , 1>1|lil] ) , 

P(S 6 ) = ^ (4E [6] + 3E [5 ,i] + 9E [4>2] + 9E [4 ,i,i] + 8E [3)3] + 10S [3 , 2l i] + 20E [3> i,i I i ] 

+31E[ 2 ,2,2] + 31E[ 2 ,2,i,i] + 67E[2,i,i,i,i] + 203E[i, 1,1, 1,1,1]) , 

P(SV) = -j (2E [7] + 5E [6 ,i] + 5E [5)2 ] + 7E [5 ,i,i] + 7E [4)3] + 15E [4)2j i] + 25E [4 ,i,i,i] 

+13E[3, 3 ,i] + 19E[ 3 ,2,2] + 27E[ 3 ,2,i,i] + 67E[3,i,i,i,i] + 59E[ 2 , 2,2,1] + 97E[ 2 , 2 , 1,1,1] 
+255E[ 2 , 1,1, 1,1,1] + 877E[i,i,i,i,i,i,i]) , 

P(5 8 ) = ^ (4E [8] + 3E [7 ,i] + 11S [6j2] + HS [6 ,i,i] + 5E [5i3] + 10E [5)2>1] + 20E [5|1>ljl] 

+ 16E[4, 4 ] + 13E[4, 3 ,i] + 38E[ 4 ,2,2] + 38E[ 4 ,2,i,i] + 82E[ 4 ,i, 1,1,1] + 21E[ 3j3)2 ] 
+33E[3, 3 ,i,i] + 43E[ 3 ,2,2,i] + 87E[ 3 ,2, 1,1,1] + 255E[ 3 , 1,1, 1,1,1] + 164E[ 2 ,2,2,2] 
+164E [2 ,2,2,i,i] + 352E [ 2,2,i,i,i,i ] + 1080E [2il|1)1)1)1)1] + 4140E [1>1)1i1 , 1)1j1j1] ) . 

(D.l) 

In (ID.2j) we give the multiplicities A4\ for various U(3) representations of the BPS 



states organized by total charge n. According to (I3.3ip they are equal to Xa(P), using 
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the elements in fID.lj) . 



n 


£7(3) representations 


1 


1 [11 


2 


2 [2] 


Q 

o 


<3 [OJ tf 1 [a, 1J 


4 


5 [4] © 2 [3, 1] © 2 [2, 2] 


5 


7 f5l © 5 [4 11 © 4 [3 21 
©1 [2,2,1] 


6 


11 [6] ©8 [5,1] ©10 [4, 2] ©2 [3, 3] 
©1 [4, 1,1] ©2 [3, 2,1] ©2 [2, 2, 2] 


7 


15 [7] © 15 [6, 1] © 17 [5, 2] © 10 [4, 3] 
©2 [5, 1, 1] © 7 [4, 2, 1] © 1 [3, 3, 1] © 4 [3, 2, 2] 


8 


22 [8] © 23 [7, 1] © 32 [6, 2] © 20 [5, 3] © 12 [4, 4] 
©5 [6, 1, 1] © 14 [5, 2, 1] © 9 [4, 3, 1] © 11 [4, 2, 2] © 2 [3, 3, 2] 



In ( ID.3j) we list combinations of R, r which are allowed at various values of A to write 
the free operators Ca,m a; r,t- The r multiplicity for any A,R is given by c(R,R,A). We 
skip A = [n] because in that case R runs over all partitions of n with multiplicity 1, and 
we will constrain ourselves to the U(2) sector A = [Ai, A2]. Note that the number of free 
operators is larger than 



A 


Operators 


A = 


[2,1] 




Cfl=[2,l],r=l 


A = 


[3,1] 


C[3,l], C[2,l,l] 


A = 


[2,2] 


C[3,l], 0[ 2 ,2], C[2,l,l] 


A = 


[4,1] 


C[4,l], ^[3,2], C[3,l,l]> C[2,2,l], ^[2,1,1,1] 


A = 


[3,2] 


C[4,l], C[3,2], C[3,l,l],r=l, C[3,l,l],r=2, C[2,2,l] , C[2,l,l,l] 


A = 


[5,1] 


C[5,l], C[4,2], C[4,l,l], C[3,2,l],r=l, C , [3,2,l],r=2, C[3,l,l,l], ^[2,2,1,1]) ^[2,1,1,1,1] 


A = 


[4,2] 


^[5,1], C[4,2],r=l> C[4,2],r=2, C*[4,1,1],t=1, C*[4,1,1],t=2, ^[3,3], ^[3,2,1] ,r=l , ^[3,2,1] ,r=2, 
^[3,2,1] ,r=3, C[3,l,l,l],r=l, C^[3,l,l,l],r=2, ^[2,2,2], C[2,2,l,l],r=l> ^[2,2,l,l],r=2, ^[2,1,1,1,1] 


A = 


[3,3] 


C[4,l,l]) ^[3,2,1]^=!, 0[3,2,l],r=2, ^[3,1,1,1] 



(D.3) 

E Integer sequences in coefficients of P 

The integer coefficients of P which we called t p contain a wealth of combinatoric data. 
The first few examples of t p can be seen in ( ID. II) . Here we explore a few known integer 
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sequences found among them. For more information about the sequences we will refer to 
[38]. 

An easily identifiable two-dimensional class of coefficients is found by considering p = 
[j n ]. For example: 

t[in] = {1,2, 5, 15, 52, 203, 877, 4140, . . .} (E.l) 

t [2n] = {2, 7, 31, 164, 999, 6841, . . .} (E.2) 

t [3n] = {2, 8, 42, 268, 1994, 16852, . . .} (E.3) 

t m = {3, 16, 111, 931, 9066, 99925, . . .}. (E.4) 

According to the formula (I3.37P we can express the generating functions for this class of 
sequences as 



t(Vj) = ex P 

where d are divisors of j. For example: 



£l( e *«-l)j, (E.5) 



t(yi) 


= exp (e Vl - 


-1) 




(E.6) 


t{yi) 


= exp ((e m 


-1) 


+ (e 2 ^ 2 - l)/2) 


(E.7) 


t(Vs) 


= exp ((e 2 ' 3 


-1) 


+ (e 3 ^ 3 - l)/3) 


(E.8) 




= exp ((e m 


-1) 


+ (e 2yi - l)/2 + (e %4 - l)/4) . 


(E.9) 



We find the sequences t^, tpn], t[3 n ], t[A n ] are identified in the [3B] as A000110, A002872, 
A002874, A141003 respectively. The whole class of such coefficients tyn] are in fact called 
"Sorting numbers" [7Tj . 

A couple of interesting special cases are worth mentioning. For j = 1 the sequence is 
the Bell or exponential numbers. It counts the total number of ways that n distinguishable 
elements can be divided into any number of indistinguishable subsets. 

For j = 2 we have a description as "number of partitions of 2n objects invariant under 
a permutation consisting of n 2-cycles" . 

It is also noteworthy that the individual factors in ( I3.37P have a combinatoric and 
group theoretic meaning. 



expi(x d -l)= J2B(d,k)^ (E.10) 

A;=0 

The numbers B(d,i) are discussed in [72] which defines them as 
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where aj(cx) is the number of i-cycles in the symmetric group element a. 
They are related to Stirling numbers of the second kind S(k,j). 



B^k)=^J2^r 1 (E.12) 



The Stirling numbers of the second kind S(k,j) [73] count the number of ways to partition 
a set of k elements into j non-empty subsets. An explicit formula is 

S(M = i£(-W'V- (E.13) 

F Fourier basis two-point function 

Here we present a derivation of (12.141) in the conventions of this paper. The Fourier basis 
operators are defined as (12 Ah . 

Oa,Ma,R,t = — p £ S R i f£' T D^{a)C1 MK m O Sia (F.l) 



a, a 



The two-point function in the trace basis can be evaluated from Wick contractions to be: 



(O t jO a , a ) = — (tr n (X^)|tr re (X 5 a)) 



1 

= £ w^^^-vw) (R2) 
= £ w^'V 1 ^). 

Here C(cr) is the number of cycles in permutation a and we used the definition (I2.15P 

Q = Y^N c{a) - n a. (F.3) 
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We can evaluate now the two-point function of operators fIF.ll) using (IF.2j) 

( Oa 2 ,Ma 2 ,R 2 ,t 2 I Oa u M Ai ,Ri ,n ) 

V^fli <^if 2 r, if 2 R 2 A 2 , r 2 r. Ri ifi Ai, n n if 2 f o\ n ifi / \ /W V p, 

( n \)2 2.^1 h h me ° ii ji mi ^iaia ^J^hjl W ^A 2 ,AfA 2 ,m 2y / U Ai,M Al ,m 1 



x ^/(/^V^) 

V d Ri d R2 sr^ o Ri R2 A 2j r 2S , ifx ifx Ai, n n if 2 f.-i^nfli /W s ) V r s 

S,a,7 

vjjgA „ i? 2 i? 2 A 2 , r 2c? ifi ifi Ai, ri n R 2( s D R 2( s D R 2 r s'XMj^- n R 1 , s 
( n |)2 2^ *a ia «2 15 ii ii mi U ki 2 Vl) U kl\ 0L ) U lhVi) j r U iih\ a ) 

3,a,"/ ^ 2 

X - D m|m 2 (^') (^A 2 ,M A2 ,m 2 ) ^Ai.AfA^mi 



Xif 2 (^)r r r q if 2 if 2 A 2 , t 2 q R 2 if 2 A 2 , n C #2 #2 A 2 , r' q if 2 if, 2 A 2 , r' 

—1 0Ai,A 2 0M Al ,AiA 2 ifi,if 2 o ijjjraa ii ji mi '1 ji ™i ° is j'2 mj 

"if 2 T , 

Xif 2 (^) 



#Ai ,A 2 ^ Ai Al ,M A , ^ifi ,if 2 ,72 
«if 2 

(F.4) 

in agreement with (I2.14p . The following steps were performed: 

• In the second line the /3 and b sums were performed using delta functions. 

• In the third line we expanded the product D^^^a^Vt) and used: 

A^(7" 1 )=^(7), (F.5) 

"if 

because f2 is a central element. We also use the following transformation property 
of the CG coefficient: 



^A 2 Va 2 ,™ 2 ~~ ^A 2 ,M A2 ,m 2 D m\ m2 (l)- ( F - 7 ) 



This is easily understood using the bracket notation, where C A Ma m is the overlap 
between V® n basis \a) and the V^ [M) ® V[ n basis \A,M A ,m): 

Ci;2 Aim =( 7 (a)|A,M A ,m> 
= (a| 7 |A,M A ,m) 

= (a|A, M A ,m')(A,M A ,m'\j\A, M A ,m) 

65 



Here it is important to note that according to our definition of 7(a) the permutation 
7 acts by the right action on the vector space \a), that is: 

( 7 (a)| = (a| 7 , | 7 (a))= 7 ^ 1 |a)- (F.9) 

In the fourth line we have performed the sums over a, a, 7 using the symmetric 
group identities: 

D?^{a)D?* 2 {a) = ^WA,*A,; 2 (F.IO) 

The Clebsch-Gordan coefficients S ^ i^i^ 3 ' T are again for coupling representations 
R\ ® i?2 <8> -R3 — > 1 and not as more usual R\® R2 Rz- They are related by: 

C -Rl ^2 i?3, T Q Rl R2; R'iT /Tp 1 o\ 

D il i 2 %3 - /J—^ h i2; is ■ 



l R 3 

Also we have used that our C% M m are defined in an orthonormal basis and thus 

E ( < ^A2,M A2 ,m 2 ) ^Ai,Ar Al) mi = ^Ai,A2^A/ Al ,M A2 ^mi,m 2 - (F-13) 



• Finally, in the last line we used the orthogonality of the Clebsch-Gordan coefficients: 

]T S t f * 3 ' ri S t f* 3 ' - = 6 TUT2 . (F.14) 

A useful reference for properties of CG coefficients is 
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